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Abstract 


I  hc  'irthngfinal  arras  dcsq;n  <»f  cxpcrimcrYis  for  process  and  equipment  optimization  is 
combined  vsiih  a  scmi-cmpirical  modelling  approach  based  on  dimensional  analysis.  This 
methodologs  is  applied  to  hi  Sim  Phosphorus  Doped  Polysilicon  Process  via  LPCVD. 
Taguchi  s  method  of  pure  empiricism  vsorks  vsell  for  the  optimization  of  four  different 
rcNponses  -  thickncNS  uniformitv  across  the  wafer,  uniformity  in  a  batch,  growih  rate,  and 
the  film  rcNistisits  for  the  proces>«  under  consideration  The  average  thickness  uniformity' 
bir  the  batch  i>  optimized  to  2  5''?  standard  deviation.  In  separate  confirmation  runs  a 
growth  rate  of  d8'n  ni;n  and  a  resistivity  of  1.25  mohm-cm  are  achieved.  However,  the 
dependence  of  Taguchi’s  method  on  the  assumption  that  the  process  is  free  of  parameter 
interactions  distorts  the  prediction  of  nominal-is-the-bcsi  signal-to-noise  ratio  for  the  film 
thickness  uniformits  across  the  wafer.  As  a  solution  to  such  problems  of  purely 
experimental  approach  a  thirteen  terms  semi-empirical  model  based  on  dimensional 
analysis  is  evolved  for  the  film  thickness  uniformity  in  the  process  under  consideration  by 
using  stepwise  regression 

The  average  percent  error  in  dimensional  analysis  model  is  about  a  third  and  the  F-test  is 
on  the  average  six  to  seven  times  better  than  that  for  the  corresponding  primitive  variables 
model.  Dimensional  analysis  model  predicts  more  reliably  because  the  confidence 
intervals  of  its  predictions  are  seen  to  be  more  than  three  times  narrower  than  that  for 
primitive  variables  model  predictions.  Also,  the  dimensional  analysis  model  overcomes  the 
problem  encountered  by  Taguchi  s  method  with  regards  to  the  importance  of  the  process 
parameter  interactions  that  the  later  tends  to  neglect. 
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Abstract 

The  onhogonal  array  design  of  experiments  for  process  and  equipment  optimization  is 
combined  with  a  semi-empirical  modelling  approach  based  on  dimensional  analysis.  This 
methodology  is  applied  to  In-Sim  Phosphorus  Doped  Polysilicon  Process  via  LPCVD. 
Taguchi's  method  of  pure  empiricism  works  well  for  the  optimization  of  four  different 
responses  -  thickness  uniformity  across  the  wafer,  uniformity  in  a  batch,  growth  rate,  and 
the  film  resistivity  for  the  process  under  consideration.  The  average  thiclmess  uniformity 
for  the  batch  is  optimized  to  1.5%  standard  deviation.  In  separate  confirmation  runs  a 

growth  rate  of  98  A°/min.  and  a  resistivity  of  1.25  mohm-cm  are  achieved.  However,  the 
dependence  of  Taguchi’s  method  on  the  assumption  that  the  process  is  free  of  parameter 
interactions  distorts  the  prediction  of  nominal-is-the-best  signal-to-noise  ratio  for  the  film 
thickness  uniformity  across  the  wafer.  As  a  solution  to  such  problems  of  purely 
experimental  approach  a  thirteen  terms  semi-empirical  model  based  on  dimensional  analysis 
is  evolved  for  the  film  thickness  uniformity  in  the  process  under  consideration  by  using 
stepwise  regression. 


The  average  percent  error  in  dimensional  analysis  model  is  about  a  third  and  the  F-tesi  is  on 
the  average  six  to  seven  times  better  than  that  for  the  corresponding  primitive  variables 
model.  Dimensional  analysis  model  predicts  more  reliably  because  the  confidence  intervals 
of  its  predictions  are  seen  to  be  more  than  three  times  narrower  than  that  for  primitive 
variables  model  predictions.  Also,  the  dimensional  analysis  model  overcomes  the  problem 
encountered  by  Taguchi's  method  with  regards  to  the  importance  of  the  process  parameter 
interactions  that  the  later  tends  to  neglect. 
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1.  INTRODUCTION 


It  is  possible  to  study  a  physical  process  or  a  product  design  in  terms 
of  either  analytical  models  or  empiricism.  Analytical  models  give 
good  extrapolation  and  physical  understanding  of  the  system  even 
without  access  to  the  modelled  process.  However,  this  mechanistic 
approach  requires  extensive  effort  for  evolving  and  solving  the 
relevant  equations,  e.g.  differential  equations  of  fluid  flow  processes. 
These  equations  may  get  extremely  complex  even  for  a  system  of 
moderate  complexity  and  still  not  capture  the  completely  relevant 
physics  of  the  system.  On  the  other  hand,  the  empirical  approach 
gives  good  accuracy  within  the  experimental  domain  but  it  usually 
requires  a  large  investment  in  terms  of  equipment  time  and  expense. 
Also,  lack  of  access  to  the  desired  experimental  environment,  which 
includes  process  equipment,  trained  operators,  and  data  acquisition 
and  processing  capabilities,  may  be  a  big  constraint  on  empirical 
approach. 


1.1  Motivation 

In  this  work  the  strengths  of  empiricism  and  analyticity  are 
combined  for  the  purpose  of  equipment  design.  The  approach  to 
equipment  design  and  optimization,  as  outlined  below,  poses  a 
special  challenge  for  the  In-situ  Phosphorus  Doped  Polysilicon 
Process  via  LPCVD  because  of  the  tradeoff  between  the  responses  of 
growth  rate,  film  thickness  uniformity,  and  resistivity  for  this 
process.  As  will  be  explained  later,  this  transport-limited  deposition 
process  is  very  sensitive  to  the  geometry  of  the  furnace  and  related 
quartzware.  Therefore,  it  is  attempted  to  incorporate  these 
geometrical  factors  of  the  LPCVD  equipment  along  with  process 
parameters  into  our  semi-empirical  model. 
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1.2  Current  Work 


A  parallel  design  ot  experiments,  based  on  the  Taguchi  Method  of 
orthogonal  arras  s  (151.  i^  used  for  designing  equipment  operating 
under  a  wide  range  of  conditions  for  the  In- situ  Phosphorus  Doped 
Polysilicon  Process  via  LPCVD.  The  equipment  design  is  realized  in 
two  ways.  Firstly,  the  data  from  a  set  of  parallel  experiments  based 
on  Taguchi  s  orthogonal  array  design  are  directly  used  to  optimize 
the  equipment.  Even  though  Taguchi's  method  allows  for  process 
and  equipment  optimization  by  using  a  small  number  of 
experiments,  this  method  makes  an  important  physical  assumption 
that  most  manufacturing  processes  are  can  be  modelled  by  main 
effects  and  some  selected  interactions.  In  order  to  overcome  this 
limitation  of  Taguchi's  method  a  second  approach  is  used  in  which 
the  empirical  knowledge  acquired  from  a  set  of  parallel  experiments 
IS  combined  with  the  partial  physical  understanding  of  the  process. 
This  approach  is  used  to  develop  a  semi-empirical  or  a  "Smart 
Responu’  Surface"  model,  based  on  dimensional  analysis,  for  the 
ratio  (known  as  S.V  ratio)  of  thickness  nonuniformity  (measured  by 
standard  deviation*  to  thickness  mean  of  the  doped  polysilicon  film. 
This  model  is  tested  to  predict  the  optimum  conditions  for  equipment 
operation.  B>  comparing  these  two  approaches  the  physical 
significance  of  the  S.\'  ratio  and.  hence,  its  utility  is  highlighted. 


2.  EXPERIMENTAL  DESIGN 


The  primary  goal  of  any  experimental  design  is  to  explore  and 
improve  the  widest  possible  domain  of  process  operation  or  product 
design  with  a  minimum  experimental  effort.  Of  the  two  methods  of 
experimentation  -  parallel  and  sequential,  orthogonal  array  parallel 
design  of  experiments  w'as  selected  to  achieve  the  goal  outlined 
before. 

The  orthogonal  array  approach  has  been  popularized  by  Japanese 
quality  expert  G.  Taguchi  [15].  His  idea  of  parameter  design  to  make 
the  product  design  and  the  manufacturing  processes  robust  to  the 
uncontrollable  variations  is  based  on  the  objective  of  setting  process 
or  product  parameters  to  values  that  minimize  the  average  squared 
deviation  of  the  response  from  its  target  value. 


2.1  Orthogonal  Array  Experimentation 

One  of  the  objectives  of  this  research  effort  was  to  design  optimal 
equipment  for  the  in-situ  phosphorus  doped  polysilicon  LPCVD 
process  by  exploring  a  large  operating  space  for  the  process 
parameters.  This  operating  space  consisted  of  the  following  eight 
factors  or  process  parameters  (see  Fig.  2.1); 

-  interwafer  distance  (w), 

-  %  hole  area  on  the  cage  (Ac), 

-  deposition  temperature  (Td), 

-  deposition  pressure  (P), 

-  inter-injector  distance  dx). 

-  total  silane  flow  (Qs), 

-  phosphine  flow  through  injector  #1  (Qpi),  and 

-  phosphine  flow  through  injector  #2  (Qp2)  • 
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Picio'i?  Distance  "w" 

Fig. 2.1:  Schematic  of  LPCVD  Process  Tube 

The  effect  of  these  process  parameters  was  studied  on  the  following 
three  responses: 

-  film  thickness  uniformity. 

-  film  growth  rate  and 

-  film  resistivity. 

The  choice  of  eight  process  parameters  was  primarily  made  on  the 
basis  of  the  past  work  done  on  the  in-situ  doped  polysilicon  process 
[4.  6.  7,  8.  9.  10].  The  fact  that  eight  parameters  were  selected  to 
be  used  in  an  eight  column  L18  orthogonal  array,  imposed  two 
limitations  on  the  design: 

-  the  parameters  had  to  be  independent  of  each  other  ;  an  J 

-  the  factor  interactions  could  not  be  studied  from  the  results. 

As  will  be  seen  in  the  final  results,  these  limitations  were  not 
severely  inhibiting.  Further,  the  semi-empirical  model  was  evolved 
to  overcome  these  limitations. 


The  settings  or  levels  for  these  process  parameters  'vere  chosen  to 
span  a  vide  but  valid  and  useful  operating  space.  The  factors  and 
their  chosen  levels  are  shown  in  Table  2.1  and  the  rationale  for  their 
selection  is  discussed  in  section  2.2.  Besides  covering  a  wide 
operating  range,  anoth.r  consideration  in  designing  the  experiments 
was  the  constraint  of  time  and  expense.  Therefore,  only  a  partial 
factorial  of  the  total  number  of  possible  combinations  of  the  factor 
levels  was  considered  for  experiments.  This  partial  factorial  design 
was  based  on  Taguchi  LI 8  orthogon.il  array*  as  shown  in  Fig.  2.2.  In 


□ 

Q 

B! 

tJp 

Tl 

Qsl 

El 

0 

IB 

nn 

n 

SI 

n: 

D 

B 

D 

2) 

■ 

D 

B 

m 

0 

0 

a 

a 

3) 

■ 

n 

m 

a 

B 

a 

a 

EQ 

r 

2 

ni 

B 

Fi 

a 

a 

5) 

B 

a 

a 

B 

a 

a 

a 

a 

6) 

■ 

a 

B 

B 

a 

a 

a 

a 

7) 

1 

3 

a 

2 

1 

B 

a 

B 

8) 

1 

3 

2. 

3 

2 

1 

[3 

I] 

9) 

1 

3 

3 

1 

3 

2 

a 

a 

10) 

1 

a 

a 

ft" 

a 

aa 

a 

■1 

a 

a 

2 

LU 

a 

a 

a 

a 

■ 

a 

a 

a 

B 

a 

ri 

a 

2  12 

a 

2 

B 

a 

a 

a 

2  12 

a 

B 

a 

a 

a 

a 

15) 

aa 

la 

n 

a 

1 

a 

a 

Bai 

2  13 

a 

3 

[2 

l3 

:  1 

ID 

2  13 

a 

n 

a 

a 

a 

a 

aa 

OB 

a 

a 

a 

la 

Fig  2  2  Taguchi  L18  Orthogonal  Array’ 


The  rf  ihc  arra>  arc  muiualU  onhogonal  m  a  combinatoric  .sense 

For  ans  pair  of  columns,  all  combinations  of  factor  lescis  occur  and  they  occur 
equal  numrscr  ol  times  This  balancing  propeny  of  the  columns  implies 
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order  to  cover  the  full  factorial  of  the  factor  level  combinations  one 
needs  to  do  2^X3^  =  4374  distinct  experiments.  However,  an 
orthogonal  array  experiment  design  provides  us  with  neccessary 
data  with  only  18  experiments,  provided  the  superposition  or 
additive  model  [17]  for  the  main  effects  of  the  parameters  holds  for 
the  process  under  scrutiny. 


Process  Parameter  Levels 


Factor 

Lovols 

1 

2 

3 

1 

Interwafer  /  \ 

d„t.n« 

7.5  mm 

1  5. mm 

— 

2 

Percent  Hole  /  *  \ 
Area  on  Cage 

05% 

3.5% 

20% 

3 

Depoaition  /j  \ 
Temperature  ' 

580“C 

600PC 

620°C 

Pressure  in  /p\ 
the  Tube 

325 

mtorr 

350 

mtorr 

375 

mtorr 

5 

Inter-Injector/ 1  x 
Distance 

^0  0  cm 

52.5cm 

55.0cm 

6 

Total  Silane  \ 

Flow  Rate 

250 

seem 

300 

seem 

350 

seem 

7 

Injector  1  Phosphine 

(Qp,) 

O&scm 

^cem 

8 

_ 

InjectorZ  Phosphine 
Flow  Rate  (Q^^) 

26,cc,r, 

07 

seem 

1.8 

seem 

Table  2.1:  Process  Parameter  Levels 


Table  2.1  shows  the  factor  levels  used  in  the  experiments.  A  limited 
range  of  two  levels.  ever>  slot  in  the  boat  (  level  1  )  and  every  other 

1  0 


slot  (  level  2  ).  was  chosen  for  interwafer  distance  "w"  because  it 
was  desired  that  the  film  properties  be  improved  without  cutting  the 
load  size  drastically.  Since  "w"  had  only  two  levels,  it  was  placed  in 
column  one  of  the  orthogonal  array  in  Fig  2.2^  The  levels  for  %  hole 
area  "Ac"  on  the  cage  were  chosen  based  on  the  past  work  done  in 
the  process  lab  at  BTU  Engineering.  According  to  this  work  the  film 
thickness  uniformity  deteriorated  rapidly  for  effective  cage  hole  area 
above  lO^c.  A  wide  range  of  0.5-20%  "Ac"  was  selected  and  divided 
geometrically  into  three  levels  shown  in  Table  2.1.  The  number  of 
levels  were  limited  to  three  because  the  quartzware  cost  for  the 
cages  is  high.  Since  changing  the  cage  between  experiments  is  time 
consuming,  the  process  parameter  "Ac"  was  placed  in  second  column 
of  the  orthogonal  array  so  that  more  experiments  could  be  completed 
with  the  same  cage  setup.  The  levels  for  deposition  temperature  "Td" 

uere  based  on  the  fact  that  the  in-situ  doped  silicon  film  via  LPCVD 

is  amorphous  when  deposited  at  temperature  below  580°C  and  the 
film  surface  roughness  deteriorates  rapidly  above  620°C.  Although 
higher  deposition  temperature  will  enhance  the  growth  rate  of  the 
film,  the  out-diffusion  of  the  dopant  molecules  become  very 
significant  at  higher  temperatures  thereby  raising  the  film 
resistivity.  Keeping  in  view  the  tradeoffs  between  film  quality, 
growth  rate  and  film  resistivity.  "Td"  levels  of  580°C,  600°C  and 
6  20^C  were  chosen.  The  levels  for  deposition  pressure  "P"  were 
determined  by  low  pressure  considerations  for  this  LPCVD  process 

and  by  the  constraint  imposed  by  the  capacity  of  the  vacuum  pump 
used.  The  vacuum  pump  failed  to  pump  pressures  below  325  mtorr 
at  the  highest  silane  flow  rate  of  350  seem,  "lx"  levels  w'ere  chosen 
based  on  the  past  work  in  the  BTU  Engineering  process  lab  for 

controlling  down-the-load  film  thickness  unformity. 

Flow  rates  for  silane  and  phosphine  through  the  two  injectors  were 
chosen  based  on  the  base  line  runs  made  in  the  past  by  the  BTU 
Engineering  staff.  The  ratio  of  the  flows  of  either  of  these  gases 


'•The  entries  "1",  "2'.  "3"  in  the  onhogonal  array  represent  the  levels  of  the 
process  parameters.  Row  numbers  1-18  represent  the  run  numbers. 


through  injector  #1  to  that  through  injector  #2  was  2. 2: 1.0.  The 
lower  flow  rates  in  injector  #2  is  explained  by  the  fact  that  the  role 
of  the  second  injector  is  only  to  compensate  for  down-the-load 
depletion  effects,  if  any.  The  total  silane  flow  through  the  two 
injectors  was  considered  as  one  parameter  rather  than  tw'o  different 
injector  flows  because  silylene,  and  not  silane,  is  the  dominant 
reactant  in  this  transport  limited  process.  Therefore,  it  is  the  total 
amount  of  silane  flowing  through  two  injectors  and  its  homogeneous 
phase  product  silylene  that  is  important  and  not  the  axial 
distribution  of  silane  in  the  furnace.  In  any  case  the  axial 
distribution  of  silane  and  the  consequent  down-the-load  depletion 
effects  are  controlled  by  "Ix".  Also,  the  coupling  of  two  silane  flows 
into  one  parameter  reduced  the  number  of  process  parameters  to 
eight  so  that  they  could  fit  in  a  LI  8  orthogonal  array.  In  contrast, 
the  phosphine  flow  rates  were  independently  controlled  so  that  their 
distribution  in  two  injectors  could  be  independently  controlled  to 
influence  across-the-wafer  film  thickness  uniformity  down  the  batch. 

2.3  Noise  Factors  and  Testing  Conditions 

Although  the  influence  of  the  process  parameters  was  studied  by 
LI 8  experimental  design,  the  influence  of  various  random  variations 
on  the  process  demanded  an  extended  approach  of  identifying 
various  noise  factors  and  studying  their  influence  on  the  responses 
under  consideration. 

Some  of  the  noise  factors  identified  were:  predeposition  anneal, 
paddle  off-center  in  the  furnace,  humidity,  room  temperature,  wafer 
thickness  variations,  impurity  atoms  in  the  lattice,  film  growth  on  the 
tube  liner,  and  measuring  device  limitations.  Except  for  the 
predeposition  anneal  and  paddle  off-center,  all  other  noise  factors 
were  not  controllable.  In  order  to  study  the  influence  of  noises  in 
the  system,  replicate  run  was  made  at  the  end  of  the  L18 
experiments.  The  average  error  variance  introduced  by  the  various 
noise  factors  (i.e.  the  difference  between  the  replicate  runs)  was 


much  smaller  than  the  main  effects  of  the  process  parameters.  It 
was  concluded  that  the  repeatablity  of  the  process  was  not  greatly 
affected  by  the  noises  in  the  system  and.  therefore,  noise  factors 
were  not  considered  in  the  subsequent  steps.  The  comparison  of  a 
replicate  with  original  run  #1  in  L18  is  shown  in  Fig. 2. 3. 
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Fig.  2.3;  Replicate  vs  Original  for  Run  #1  in  LI  8 
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2.4  Data  Analysis  Plan 

The  noise  factors  and  random  fluctuations  in  the  process  parameters 
vary  the  quality  characterstic  of  a  product  from  unit  to  unit  within  a 
batch.  The  average  of  squared  deviation  is  called  the  quadratic 
quality  loss.  Such  quadratic  quality  loss  is  expressed  in  dollar  terms 
by  multiplying  the  quality  loss  with  a  loss  coefficient. 


2.4.1  Cost  Function 


In  order  to  optimize  the  thickness  uniformity,  growth  rate  and 
resistivity  of  the  polysilicon  film,  a  total  cost  function  was  derived 
that  incorporated  the  quality  loss  due  to  inhomogeneous  films,  low 
growth  rate  and  high  resistivity.  However,  lack  of  data  for  quality 
cost  coefficients  (defined  below’)  prevented  us  from  calibrating  the 
total  cost  function.  Since  these  responses  are  in  tradeoff  with  each 
other,  their  individual  optimization  was  also  done  in  order  to 
compare  their  results  and  highlight  the  tradeoff. 


Accordingly,  the  following  loss  functions  were  derived  for  the 
concerned  responses  (see  Appendix  B  for  derivation); 


1)  Thickness  Uniformity  Cost  :  Film  thickness  mean  and 

standard  deviation  "OT,i"  was  observed  for  each  ith  wafer  in  every 
batch.  Based  on  these  observations  two  different  adjusted  quality 
loss  [14,  15,  17]  functions  were  considered:  across-the-wafer  quality 
loss  and  overall  loss  in  a  batch.  Considering  nominal-the-best 
quadratic  quality  loss  function,  the  total  quality  loss  due  to  across- 
the-wafer  thickness  nonuniformity  after  the  mean  is  adjusted  to  the 
target  value  is  just  the  sum  of  losses  due  to  individual  wafers  and  is 
given  by. 


QLti  —  kjj 
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(2.1) 


where  kji  is  the  film  thickness  quality  loss  coefficient  for  a  wafer,  nw 
is  the  number  of  wafers  in  the  batch  and  is  the  target  thickness. 
The  corresponding  signal-to-noise  ratio  SNxi,  as  defined  by  G. 
Taguchi,  is  given  by. 
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Note  that  the  constants  kji  and  do  not  appear  in  equation  2.2 
because  thay  do  not  vary  from  one  run  to  another  and,  therefore, 
they  play  no  role  in  the  optimization.  The  logarithmic  transformation 
expresses  the  quality  loss  in  a  decibel  scale. 


Similarly,  total  nominal-the-best  quadratic  quality  loss  due  to 
thickness  variance  across-the-wafer  and  across-the-load  is  given  by, 


where  kT2  is  the  film  thickness  quality  loss  coefficient  for  a  batch,  M-j. 
is  the  mean  thickness  for  the  batch  and  oj  is  the  overall  standard 
deviation  for  the  batch  (see  Appendix  B).  The  signal-to-noise  ratio 
SNj:  for  equation  2.3  is  given  by. 


SN,,=  -10log,„(4’ 


(2.4) 


.4s  before,  constants  kT2  and  <io  not  appear  in  equation  2.4  , 

2)  Growth  Rate  Cost:  An  enhancement  in  growth  rate  of  the  film 
leads  to  shorter  deposition  time  for  achieving  target  thickness.  The 
cost  function,  as  discussed  in  Appendix  B,  is  given  by. 
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and  the  corresponding  decibel  scale  signal-to-noise  ratio  SNg  is  given 
bv. 


SN<3=  -lOlog., 
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where  kc  is  the  growth  rate  cost  coefficient,  and  to  is  the  deposition 
time.  The  constant  terms  in  equation  2.5  do  not  appear  in  equation 
2.6  for  the  reason  discussed  before. 


It  is  important  to  note  that  the  cost  function  in  equation  2.5  is  not  a 
cost  due  to  quality  loss  because  the  load  is  never  unacceptable  even 
if  the  growth  rate  falls  indefinitely.  To  counter  the  fall  in  growth 
rate  the  adjustment  factor  of  deposition  time  is  increased,  thereby 
incurring  increased  fixed  costs  for  labor,  floor  space,  and  equipment. 
The  cost  Cg  consists  of  the  labor  cost  Cl,  the  floor  space  cost  Cf,  the 
equipment  cost  Ce- 


3)  Resistivity  Cost:  The  target  was  to  minimize  the  resistivity  of  the 
film  by  controlling  the  process  parameters.  Resistivity  is  a 
continuous  and  non-negative  quality  characterstic  with  no 
adjustment  factor  in  the  deposition  process.  Therefore,  we  simply 
minimized  the  quadratic  quality  loss  without  adjustment  [14],  given 
by. 


QLp  = 


(2.7) 


w'here  j  is  the  mean  resistivity  across  the  ith  wafer  in  a  batch,  and 
kp  is  the  resistivity  quality  cost  coefficient.  Minimizing  QLp  in 

equation  2.7  is  equivalent  to  maximizing  smaller-the-better  signal- 
to-noise  ratio  SXp  given  by. 


SN  =  -101o2. 
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As  before,  we  ignored  the  constant  kp  and  expressed  the  quality  loss 
in  decibel  scale. 


4)  Total  Quality  Cost:  In  order  to  find  an  optimal  point  for  the 
process  parameters  with  the  constraint  of  the  tradeoff  between  the 
three  responses  considered,  i.e.  the  growth  rate,  film  thickness 
uniformity,  and  the  resistivity,  a  total  quality  loss  function  was 
considered.  It  was  the  equally  weighted  sum  of  the  costs  in 
equations  2.3,  2.5  and  2.7  and  is  given  by, 

QL„v«„a  =QLt2+Cc+QLp  (2.9) 

This  total  quality  cost  is  expressed  in  a  decibel  scale  to  give  an 
overall  signal-to-noise  ratio, 

SN„„.,  =  -lCllog,jQL,,+Q+Qg 


(2.10) 


The  target  was  to  minimize  the  total  cost  in  equation  2.9  by 
maximizing  the  SNoveraii  in  equation  2.10  . 


The  cost  coefficients,  as  determined  in  Appendix  B,  are  given  by. 
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(2.11(b)) 


(2.11(c)) 


where  Ao  is  the  S  cost  incurred  in  discarding  the  entire  batch  when 
the  batch  non-uniformity  exceeds  A:  angstroms  (A°)  for  a  target 

thickness  of  ^,t  A°,  or  resistivity  overshoots  the  functional  limit  of 
.cm  on  a  target  resistivity  of  r2_cm  (Note  that  ideally  this 
target  is  zero).  The  terms  cl,  cf,  ce  represent  the  fixed  cost  rate 
(S/min)  for  labor,  floor,  and  equipment  respectively.  Also  the  total 

predeposition  and  postdeposition  steps  time  t]  is  fixed  for  all  runs. 

« 

Time  represents  the  time  it  takes  for  the  deposition  step  to  achieve 
target  thickness. 

2.4.2  ANOVA  Analysis 

After  doing  the  cost  function  analysis  of  the  observed  data,  the  next 
step  was  to  estimate  the  main  effects  of  the  process  parameters  and 
to  perform  analysis  of  variance  (ANOVA)  [16].  Analysis  of  variance 
of  the  data  gave  the  relative  effect  of  the  process  parameters  on  the 
response  considered.  From  this  the  most  significant  main  effects 
were  determined.  ANOVA  was  performed  for  all  the  different  cases 
of  SN  ratios  discussed  above. 


3.  **Smart  Response  Surface*^ 

MODELLING 


PROCESS 


Data  have  no  meaning  in  themselves;  but  they  assume  meaning  in 
relation  to  a  conceptual  model  of  the  phenomenon  under  scrutiny.  In 
this  chapter  it  is  attempted  to  utilize  the  partial  understanding  of  the 
in-situ  doped  polysilicon  process  and  combine  it  with  the  empirical 
data  to  evolve  a  semi-empirical  or  a  "Smart  Response  Surface"  model 
based  on  dimensional  analysis.  The  underlying  motivation  for  the 
semi-empirical  approach  is  to  do  without  the  intricacies  of 
mechanistic  modelling,  such  as  the  solving  of  three  dimensional 
differential  equations  for  mass  transport  and  reaction  kinetics, 
without  sacrificing  the  causal  accuracy  of  the  predictions  of  the 
resultant  model.  This  causal  accuracy  is  achieved  by  combining  the 
process  variables  into  dimensionless  groups  which  are  physically 
more  intuitive  than  the  primitive  process  variables. 


3.1  Process  Background 

Low  pressure  chemical  vapor  deposition  (LPCVD)  in  hot  walled 
reactors  is  rapidly  replacing  the  conventional  atmospheric  pressure 
CVD  (APCVD)  method  of  depositing  films  on  the  w'afers  for  various 
applications  in  microelectronics  technology.  The  LPCVD  method  for 
producing  thin  insulators  and  semiconductor  films  gives  reduced 
processing  cost  and  superior  product  quality  in  terms  of  fewer 
defects,  better  thickness  uniformity,  and  substantially  improved  step 
coverage  than  APCVD  films.  The  close  spaced,  stand-up 
configuration  for  device  wafers  permissible  in  LPCVD  reactors 
(Fig.3.1)  leads  to  greatly  increased  throughput'  despite  the  reduced 
film  deposition  rates. 


'load  size  divided  by  deposition  cycle  time 


Liner  Wafers  Process  Quartz 


The  LPCVD  technique  for  depositing  polycrystalline  silicon  films  by 
the  pyrolysis  of  silane  (SiH4)  yields  very  uniform  films  because  at 
reduced  pressure  diffusion  is  rapid  enough  so  that  the  chemical 
reactions  at  the  wafer  surface,  and  not  the  gas  phase  mass-transfer 
process,  are  the  rate  controlling  factorsf2,  11].  However,  the 
application  of  LPCVD  to  the  in-situ  phosphorus  doped  polycrystalline 
process  degrades  the  polysilicon  LPCVD  process  because  of  the 
phosphine  gas  (PH3)  introduced  into  the  furnace  along  with  silane,  [4, 
5].  This  logical  extension  of  the  intrinsic  polysilicon  process  is 
desirable  because  it  avoids  the  complexity  of  more  conventional 
doping  techniques  such  as  post-deposition  implantation  of  dopant,  or 
indiffusion  of  dopant  from  a  dopant  rich  layer  applied  on  the  silicon 
surface,  and  offers  lower  film  resistivity  because  of  more  uniform 
incorporation  of  phosphorus.  But  the  introduction  of  phosphine  into 
the  system  severly  suppresses  the  growth  rate,  leads  to  wiihin-wafer 
and  wafer-to-wafer  film  thickness  non-uniformities,  and  poses  a 
problematic  tradeoff  between  growth  rate,  thickness  uniformity,  and 
resistivity  [8,  9]. 

The  most  widely  accepted  explanation  for  the  problems  of  growth 
rate  suppression  and  film  inhomogeneties  is  that  phosphine 
passivates  the  silicon  surface  by  altering  the  silicon  growth  kinetics 


[4.  5).  A  lone  pair  of  electrons  on  the  phosphine  molecule  helps  it  to 
compete  successfully  with  silane  (a  fully  coordinated  molecule)  for 
the  active  surface  sites  and  form  a  layer  of  nondissociative 
pentavalent  phosphorus  bridging  with  adjoining  silicon  atoms.  This 
exclusion  of  silane  from  the  heterogeneous  chemistry  leads  to 
tremendously  lower  growth  rates.  It  has  been  shown  that  the 
adsorption  of  phosphine  proceeds  about  40  times  faster  than  that  of 
silane  [4], 

The  second  problem  of  film  non-uniformities  renders  the  films  unfit 

for  use  in  the  VLSI  applications.  The  films  deposited  in  this  process 
have  highest  growth  rate  along  the  edge  of  the  wafer  and  reduce  to  a 
minimum  at  the  wafer  center  (Fig3.2).  This  growth  rate  gradient  is 

directly  related  to  the  concentration  gradient  of  the  active  species  in 
the  interwafer  region.  Also,  the  active  species  participating  in  the 
film  deposition  is  no  longer  silane  (it  is  suppressed  by  phosphine) 

but  silylene  (SiH:).  a  homogeneous  species  resulting  from  the 

pyrolysis  of  silane 


SiH^Cg)  < . >  SiH:(g)  -i-  H:(g) 

Silylene  has  an  activation  energy  comparable  to  that  of  phosphine 
(both  have  a  lone  pair  of  electrons)  and,  therefore,  can  successfully 
insert  into  the  Si-H  bonds  on  the  silicon  surface.  The  suppression  of 
the  heterogeneous  chemistry  of  silane  (a  majority  homogeneous 
species)  by  phosphine  and  the  fractional  contribution  of  silylene  (a 
minority  homogeneous  species)  to  the  growth  rate  indicates  that  this 
process,  in  contrast  with  the  intrinsic  polysilicon  LPCVD,  is  mass 
transport  limited  and  not  surface  kinetic  controlled  [11].  Thus,  the 
depletion  of  the  active  silylene  species  in  the  interwafer  region 
results  in  the  growth  rate  profiles  shown  in  Fig3.2.  This  depletion 
effect,  or  bull's  eye  effect  -  as  it  is  called,  can  also  be  looked  upon  as 
the  silylene  contribution  from  the  annular  region  depositing  on  the 
wafer  edges.  To  exclude  the  annular  silylene  from  interwafer  region 
and,  therefore,  improve  the  growth  rate  uniformity,  it  is 
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recommended  that  the  wafer  load  be  placed  in  a  cage  with  a  surface 
area  large  enough  to  consume  the  annular  silylene. 


ootated  8 
lilted  30 

6.59 


Fig.  3.2;  Typical  Growth  Rate  Profile  in  In-Situ  Doped 
Polysilicon  Process 

3.2  Dimensional  Analysis  Model  Formulation 

A  model  is  developed  for  the  tradeoff  between  the  growth  rate  mean 
F  and  nonuniformity  (measured  by  standard  deviation  o)  by 
considering  the  following  "signal-to-noise"  (SN)  ratio  as  the  desired 
response  for  each  individual  wafer  in  a  batch 


SN=  -lOlogio 


(3.1) 


The  Slated  goal  is  to  optimize  the  process  under  consideration  i.e. 
increase  SN  ratio  by  means  of  a  parallel  design  of  experiments  (DOE)> 
and  compare  the  DOE  predictions  with  that  of  the  SN  model.  In  the 
process  of  doing  so.  the  strength  of  dimensional  analysis  as  a  Smart 
Response  Surface  methodology  will  be  demonstrated. 

Consider  the  "reaction  cell"  shown  in  Fig. 3.?.  Of  all  the  geometrical 
parameters  shown,  only  interwafer  distance  "w"  and  the  cage  hoie 
radius  "rv."  are  the  geometrical  variables.  From  here  on  all  other 
constant  length  parameters  are  referred  to  as  "L,"  . 
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Fig. 3  3  Schematic  of  the  "Reaction  Cell" 


The  following  process  parameters  (with  their  dimensions)  are 
identified  and  considered: 

(Process  Variables) 

-  Interwafer  distance  (w) 

-  Hole  radius  on  the  cage  (rj,) 

-  Deposition  temperature  (Tj) 

-  Deposition  pressure  (P) 

-  < 


(LI. 

[L]. 

[e], 

(M/LT^l. 


'  See  Chapter  2  for  a  diecussion  on  S.N  ratios 


-  Inier-injecior  distance  (I^)  [L]. 

-  Total  silane  flow,  (Qs»  [M/T]. 

-  Phosphine  flow  through  injector  #1  (Qpi)  [M/T], 

-  Phosphine  flow  through  injector  #2  (Qp:)  [M/T], 

-  Multicomponent  Gas  Diffusivity  (Dij)  [L-/T]. 

(Process  Constants) 

-  Constant  geometrical  parameters  (Lj)  [L], 

-  Ideal  Gas  Constant  (Rq)  [ML^/T-B], 

-  Actuation  Energies  (Ej)  [ML2/T2]. 

-  Molecular  weights  of  various  species  (Mj)  [M], 


-  Film  surface  densities  of  various  species  (Pi)[M/L'^]. 

(Process  Response) 

-  Response  Parameter  (SN)  [*0-]. 

The  diffusis  ity  D.,  is  defined  by  Tj.  P.  M,.  and  Rq  and.  therefore,  can 
be  neglected.  Similarly.  P;  is  neglected  because  it  is  defined  by  M, 
and  L..  .At  this  stage  we  have  12  process  parameters  (8  variables 
plu''  2  constants). 

llasing  listed  the  significant  process  parameters,  the  Pi  Theorem  of 
Dimensional  Analysis  is  utilized  to  develop  a  set  of  dimensionless 
\ariables  called  the  Pi  Groups.  These  Pi  groups  must  satisfy  the 
condition  of  dimensional  homogenciiy  i.e.  all  the  Pi  groups  in  an 
additive  model  must  have  the  same  dimensions.  Since  the  response 
SN  is  dimensionless,  all  other  Pi  groups  must  also  be  dimensionless 
for  any  meaningful  solution  to  the  process.  From  the  above 
parameter  list,  one  parameter  (variable  or  constant)  for  each 
dimension  (length  L.  time  T.  mass  M.  and  temperature  6)  is  selected, 
making  sure  that  the  parameters  chosen,  called  the  primary 
parameters  from  hereon,  are  dimensionally  independent  i.e.  they 
cannot  be  power  grouped”  into  a  dimensionless  Pi  group.  The 
pritTiars  parameters  selected  are; 


-  Constant  geometrical  parameter  (Lj)  for  length  dimension. 


-  Molecular  weight  (Mi)  for  mass  dimension. 

-  (  Ro)/( Ej)  for  temperature  dimension. 

-  Pressure  (P)  for  time  dimension. 

It  can  be  shown  that  there  is  no  power  grouping  of  these  primary 
parameters  which  is  dimensionless.  Even  though  we  can  cancel  mass, 
length,  and  time  dimension,  no  power  group  involving  these  primary 
parameters  rids  itself  of  the  temperature  dimension.  Also,  having 
selected  5  of  the  above  mentioned  12  process  parameters  as  the 
primary  parameters.  Pi  theorem  states  that  the  number  of  Pi  groups 
should  be  no  more  than  7. 

Now  comes  the  most  critical  step  of  using  these  primary  parameters 
to  nondimensionalize  other  parameters,  called  secondary  parameters 
from  hereon.  This  is  where  the  secondary  parameters  are  grouped 
with  primary  parameters  on  the  basis  of  some  understanding  about 
the  physical  reality  of  the  process  under  consideration.  Only 
secondary  process  variables  need  to  be  nondimensionalized  because 
the  process  constants  merge  with  a  common  constant  term  Cq  in  an 
additive  model  for  the  process  response. 

Geometrical  variables  are  nondimensionalized  by  L,  to  give  the 
following  Pi  groups. 
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Deposition  temperature  will  interact  with  the  activation  energies  to 
determine  the  rate  of  silylene  production  in  the  gas  phase  and  also 
the  adsorption  of  various  species  onto  the  polysilicon  film.  This  gives 
us  a  Pi  group. 
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The  variation  in  depostion  temperature  is  expected  to  influence  the 
concentrations  in  gas  phase,  given  by  (P/RoTd),  but  this  influence  is 
found  to  be  very  small  because  of  a  very  limited  range  of  pressure 
and  temperature  levels  allowed  for  this  process.  The  flow  rates  are 
nondimensionalized  by  pressure  to  give  the  following  Pi  groups, 
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For  the  process  under  consideration,  the  fraction  of  phosphine  to 
silane,  and  not  the  phosphine  flow,  is  the  controlling  parameter. 
Therefore,  equation  3.6  is  divided  by  equations  3.7  and  3.8  to  obtain 
physically  more  significant  Pi  groups. 


(3.9) 

(3.10) 


Based  on  these  Pi  groups  the  functional  relation  can  be  expressed  as 

SN  =f  (n,.n;,n3,n4./75,r?6./T,')  on. 

This  functional  relation  indicates  two  advantages  of  dimensional 
analysis.  Firstly,  the  number  of  process  parameters  is  reduced  from 
12  (or  8  process  variables)  to  seven  Pi  groups.  Therefore,  for  the 
future  usage  of  the  model  described  below  we  need  to  do  fewer 
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numbers  of  experiments  because  the  number  of  variables  have 

reduced.  Secondly,  the  Pi  groups  above  are  physically  more 
significant.  This  infusion  of  physical  knowledge  is  enunciated  further 
below. 

An  appropriate  functional  relation  was  reached  by  starting  with  a 

first  order  polynomial  of  all  Pi  groups.  All  these  first  order  main 

effects  were  retained  in  the  subsequent  steps.  To  this  model  various 
second  order  terms  were  added  and/or  deleted  by  using  stepwise 

regression  with  T-tests'  as  the  deciding  criterion.  Then  the  partial 
understanding  of  the  physics  of  the  process  was  used  to  include  some 
physically  sensible  interactions.  Finally,  for  the  purposes  of 
regression  evaluation  of  an  additive  model  for  SN  ratio,  the  following 
13  term  equation  was  used. 


SN=  C(,  +C,(w)  +C.(rJ  +  C3(TJ  +  Q(IJ  + 


(3.12) 

Note  that  the  constant  parameters  shown  in  equation  3.2  to  3.11 

merge  with  the  calibration  coefficients  Co,  Cj .  and  Ci:  in  equation 

3.12.  Descriptive  use  of  this  regression  model  emphasized  precise 
estimation  of  regression  coefficients  by  using  the  data  from  L18 
orthogonal  array  experiments,  and  predictive  use  focused  on  the 
minimization  of  the  prediction  errors. 


'  The  T-lcst  says  that  a  regression  coefficient  Cj  is  significant  i.e.  not  zero  if  : 

111  >  Ul-k/Z;  d)  where . 

t  =  (Calculated  Value  of  C,)/(mean  square  error  in  that  Value  of  Cj); 
k  is  the  confidence  level;  and 
d  is  the  number  of  degrees  of  freedom  for  Cj 
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Various  interactions  included  in  equation  3.12  require  some 
discussion.  Deposition  temperature  interacts  strongly  with  the  major 
ph  osphine  fraction  flowing  through  injector  #1  because  it  controls 
the  relative  adsorption  and  desorption  of  phosphine  onto  the  wafer 
surface.  The  interaction  of  temperature  with  phosphine  fraction 
through  injector  #2  is  found  to  be  not  very  significant. 

The  interaction  of  Rh  with  (Qs/Qpi)  and  (Qs/Qpz)  was  expected  to  be 
significant  because  it  determines  the  amount  of  flow  through  the 
cage.  However,  this  interaction  w'as  neglected  because  stepwise 
regression  results  indicated  otherwise. 

.As  one  moves  away  from  the  load  end  of  the  tube  the  axial  depletion 
effects  become  significant.  Therefore,  Ix  interacts  strongly  with 
(Qs/Qpi)-  However,  near  the  load  end  there  are  no  significant  axial 
depletion  effects  and,  therefore,  the  interaction  of  Ix  with  (Qs/Qpi) 
through  injector  #1  (fixed  in  position)  is  not  expected  to  be 
significant.  The  presence  of  this  interaction  in  the  above  model  for 
the  load  end  wafers  may  distort  the  corresponding  predictions  to  a 
certain  extent.  The  phosphine  fraction  flowing  through  injector  #2 
plays  an  important  role  with  regard  to  the  position  of  the  adjustable 
injector  #2  because  the  position  of  injector  #2  determines  if  the 
phosphine  flowing  through  it  will  remain  trapped  in  a  homogeneous 
phase  complex  called  monosilylephosphine  (SiH3PH2).  hence,  the 
interaction  of  Ix  with  (Qs/Qpz)-  This  complex  is  formed  by  the  gas 
phase  reaction  of  silylene  and  phosphine  [4.  5]  and  is  unstable  above 
5  5  0'’C  in  the  presence  of  surfaces.  Therefore,  if  injector  #2  injects 
the  gases  at  a  point,  say.  beyond  the  heated  load  zone,  i.e  the  region 
between  load  end  of  the  tube  and  wafer  #1,  then  this  complex  will 
tend  to  trap  more  of  phosphine  from  injector  #2  because  there  is 
lesser  surface  area  and  lower  temperature  in  this  region. 
Entrapment  of  more  phosphine  should  lead  to  better  uniformities 
across  the  wafer  and  higher  growth  rates  (especially  near  the  load 
end)  in  this  transport  limited  process  (see  the  behavior  of  SN’ti  and 
S.N’g  in  Chapter  5).  This  and  other  observations  mentioned  above 
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were  confirmed  by  the  behavior  of  equation  3.12  to  the  experimental 
data  (see  Chapter  5). 

The  dimensional  analysis  model  was  compared  to  the  corresponding 
primitive  variables  model  which  was  evolved  by  using  the  step 
regression  method  discussed  before.  All  the  first  order  main  effects 
and  the  interactions  corresponding  to  the  ones  in  equation  3.12  were 
considered  in  this  primitive  variables  model.  Therefore,  a  13  terms 
model  for  a  design  involving  8  primitive  parameters  was  determined 
as. 


SN  =  Cq  +Ci(w)  +C2(rf,) +C3(TJ  +  CjP)  + 

Q(I>CJQJ  +  qOp,)  +C8(Qp,KC9(lx)  + 

C,o(T,)(Qj  +  C„(IJQ,)  +  C,,(I,)(Qp,+Qp,) 

(3.13) 

Note  that  the  interwafer  distance  "w"  is  regressed  only  to  a  linear 
term  because  it  has  only  two  levels  in  the  LI 8  orthogonal  array. 
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4.  EXECUTION  OF  EXPERIMENTS 


Based  on  the  experimental  design  discussed  in  the  chapter  2,  the 
experiments  were  performed  in  the  process  lab  at  BTU  Engineering. 
Before  exposing  the  wafers  to  the  in-situ  Phosphorus  Doped 
Polysilicon  Process,  each  wafer  was  treated  with  a  wet  oxidation 
process  to  grow  a  1000  A°  oxide  film  on  <100>  type  silicon  substrate. 


4.1  Equipment  Description  and  Procedure 

The  film  deposition  was  done  in  a  commercial  BTU  Engineering/Bruce 
Systems  735 1C  horizontal,  hot-wall  diffusion  furnace  with  three  zone 
temperature  profile  control  maintaining  the  desired  temperature  in  a 
32  inch  long  load  zone(see  Fig.  2.1).  The  furnace  tube  was  lined  with 
a  quartz  liner.  The  total  production  load  of  114  wafers,  placed  in  a 
32  inch  flat  temperature  profile  load  zone  on  a  paddle,  consisted  of 
three  quartz  boats  with  38,  4  inch,  575  micron  thick  wafers  each.  Of 
these  114  wafers.  12  were  test  wafers  and  then  rest  were  dummy 
wafers. 


Fig. 4.1:  Test  Wafers  Position  in  the  Load 


The  position  of  the  test  wafers  is  shown  in  Fig. 4.1.  The  interwafer 
distance  "w",  as  discussed  in  the  last  chapter,  had  two  levels  -  every 
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slot  in  the  boat  (7.5  mm  apart)  or  every  other  slot  in  the  boat  (15 
mm  apart).  The  deposition  pressure  was  held  at  a  constant  value  by 
regulating  the  pump  cross-section  with  a  butterfly  valve  controlled 
from  a  feedback  control  circuit.  The  inter-injector  distance  "lx"  was 
changed  after  each  run  by  moving  the  adjustable  injector  #2.  Three 
sets  of  quartz  cages  with  different  hole  areas  were  used  and  their 
ends  w'ere  covered  shut  with  quartz  baffles.  Two  silane  mass  flow 
controllers  with  300  seem  capacity  each  and  two  phosphine  mass 
flow  controllers  with  a  capacity  of  10  seem  (for  injector  #1)  and  6.2 
seem  (for  injector  #2)  were  used. 


Tim»  (min) 

Fig. 4. 2:  In-situ  Phosphorus  Doped  LPCVD  Process  Recipe 

The  7351  microprocessor  recipe  used  to  run  the  furnace  is  shown  in 
Fig. 4. 2.  The  process  parameters  of  step  #7  in  the  recipe  above  were 
varied  from  one  run  to  another  but  the  same  recipe  was  followed 
with  a  deposition  run  of  two  hours.  The  eighteen  experiments 
outlined  in  the  LI 8  orthogonal  array  in  Fig. 2. 2  were  performed  in  a 
random  order  so  that  the  statistical  bias  was  minimized. 


The  processed  wafers  were  then  annealed  at  1000°C  for  30  minutes 
in  a  slightly  oxidizing  atmosphere  of  95  vol  %  Sj  and  5  vol  %  Oi  in 
order  to  prevent  the  outdiffusion  of  the  dopant  phosphorus.  The  thin 
SiO;  film  formed  during  annealing  was  removed  by  7:1  buffer  oxide 
etch  before  data  acquisition. 

Test  wafers  #2.  #7,  and  #11  were  then  measured  for  film  thickness 
with  a  Prometrix  SM2000  SpectraMap  system.  Film  thickness  mean 
and  variance  for  a  wafer  was  obtained  from  49  readings  at  points 
shown  in  Fig. 4. 3  The  required  refractive  index  adjustment  factor 
(which  varies  with  the  dopant  concentration  and  the  deposition 
temperature)  was  obtained  for  each  run  from  correlation  with  the 
Dektak  surface  profilometer  thickness  readings.  This  adjustment 
factor  was  consistently  found  to  almost  1.0.  Therefore,  the  influence 
of  the  variation  of  refractive  index  on  thickness  readings  from  the 
Prometrix  was  concluded  to  be  small  and.  hence,  neglected  for  the 
purpose  of  the  decibel  scale  SN  ratios.  The  film  resistivity  was 
obtained  from  the  Prometrix  automatic  four  point  probe  readings  for 
sheet  resistance. 


Fig. 4. 3:  Typical  Thickness  Profiles  with  the  Points  of  Measurement 
by  Prometrix  SM2000  Spectramap  System  on  a  4  inch  Wafer 


5.  RESULTS  AND  DISCUSSION 


The  observed  data  is  tabulated  in  Appendix  A  in  Tables  A. 1(1)  to 
A.  1(1 8).  The  measured  data  was  subjected  to  the  analysis  outlined 
in  chapter  2.  The  process  parameters  in  the  LI 8  design  were 
optimized  according  to  signal-to-noise  (SN)  ratios  outlined  in 
equation  2.2,  2.4,  2.6  and  2.8: 


The  values  for  these  SN  ratios  are  shown  in  Table  5.1  and  the 
corresponding  averaged  main  effects  for  the  process  parameters  are 
shown  in  Table  5.2.  These  main  effects  for  a  particular  factor  level 
setting  were  determined  by  averaging  the  SN  ratio  for  the 
experiments  in  which  that  particular  factor  level  was  used.  For 
example,  the  SN  ratios  for  experiments  1  to  9  were  averaged  to 
obtain  the  main  effect  (in  decibels)  of  "w"  at  level  1  on  the 
corresponding  response.  The  results  in  Table  5.1  and  5.2  were  used 
to  decide  upon  the  optimum  levels  of  the  process  parameters  based 
on  the  criterion  of  maximizing  the  corresponding  SN  ratio. 
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1  SN 
\  ratios 

Run 

SNti 

SNt2 

SNg 

SNp 

1. 

46.992 

18.118 

35.990 

-22.126 

2. 

31.712 

26.545 

34.156 

-18.190 

3. 

25.493 

18.124 

34.090 

-18.848 

4. 

24.647 

20.569 

30.132 

-19.388 

5. 

40.843 

20.215 

38.662 

-31.477 

6. 

26.250 

23.123 

35.407 

-20.385 

7. 

25.455 

11.285 

33.412 

-17.251 

8. 

21.246 

17.299 

33.747 

-17.717 

9. 

31.383 

15.597 

37.912 

-24.884 

10. 

37.344 

22.045 

36.040 

-18.996 

11. 

27.794 

22.159 

34.459 

-17.358 

12. 

29.745 

11.358 

36.112 

-19.363 

13. 

24.804 

21.563 

33.202 

-17.339 

14. 

25.936 

10.210 

36.631 

-17.110 

15. 

27.719 

23.683 

38.060 

-21.258 

16. 

31.473 

24.642 

37.025 

-18.530 

17. 

22.373 

12.458 

35.485 

-16.931 

18. 

21.762 

18.990 

36.921 

-17.062 

Table  5.1:  SN  ratios  for  Orthogonal  Array  Runs 


5.1  Film  Thickness  Profile  and  Growth  Rate  Analysis 

The  main  effects  of  the  process  parameters  on  the  overall  batch 
variance  (i.e.  SNt:)  and  across-the-wafer  variance  (i.e.  SNti),  and 
growth  rate  (i.e.  SXg)  are  plotted  in  the  marginal  plots  in  Fig. 5. 3  and 
5.4.  The  optimum  factor  level  settings,  corresponding  to  the  highest 
SX's,  are  shown  in  Fig. 5. 3  and  5.4.  Tables  A. 2,  A. 3, 
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1  SN  Ratios 

SN„ 

SN,, 

SNq 

SNp 

Factor 

L-  eiS 

W 

1 

30.447 

18.986 

34.834 

-21.137 

2 

27.661 

18.567 

35.993 

-18.218 

B 

1 

33.180 

19.725 

35.141 

-19.147 

2 

19.894 

mmm 

mmm 

3 

25.616 

16.714 

-18.726 

1 

31.786 

19.706 

34.300 

-18.939 

2 

H^SISB 

mmm 

-19.794 

HSI^H 

P 

1 

30.151 

mmm 

-20.327 

2 

18.328 

wmm 

-20.114 

3 

27.957 

19.241 

35.490 

-18.592 

B 

1 

29.134 

17.317 

^ESfflBI 

-18.549 

2 

34.872 

mmm 

? 

35.898 

-21.415 

Qs 

1 

28.568 

17.320 

34.991 

-18.973 

28.797 

18.993 

35.299 

-19.272 

3 

29.797 

20.020 

35.951 

-20.787 

^pi 

1 

34.395 

16.690 

37.056 

-22.248 

2 

28.476 

19.859 

35.427 

-18.838 

3 

24.291 

19.784 

33.759 

-17.946 

Qp2 

1 

32.651 

20.058 

36.570 

-21.433 

2 

28.903 

2'>  272 

35.360 

-19.448 

3 

25,608 

14  003 

34.310 

-18.151 

Table  5.2;  Main  Effects  of  the  Process  Parameters 


and  A. 4  in  Appendix  A  show  the  ANOVA  results  for  SNti,SNt2.  and 
SN’g  respectively.  ANOVA  results  indicate  that  phoohine  flow  rates 
are  the  most  significant  parameters  influencing  the  film  thickness 
profile  and  the  growth  rate.  Batch  variance  is  also  strongly 
influenced  by  cage  hole  area,  inter-injector  distance,  and  silane  flow 
rate.  Variance  across  the  wafer  depends  on  cage  hole  area, 
deposition  temperature,  and  interwafer  distance  as  well.  Growth 
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rate  is  influenced  significantly  by  deposition  temperature,  and 
interwafer  distance  as  well. 


Fig.  5.3:  Marginal  Plots  for  SNn 

It  was  observed  that  the  film  thickness  uniformity  degrades  (i.e. 
SNti  in  Fig.  5.3  and  SXt:  in  Fig.  5.4  falls)  with  rising  interwafer 
distance  "w".  This  is  explained  by  the  fact  that  wider  "w"  means  that 
more  silylene  (SiH;^  per  wafer  is  coming  into  the  interwafer  region 
from  the  tube-cage  annulus.  This  "excess"  SiH2  deposits  near  the 
wafer  edges  to  cause  higher  nonuniformity.  This  also  leads  to 
slightly  increased  average  growth  rate  and.  therefore,  the  SN’g  rises 
with  increasing  "w"  in  Fig. 5. 5. 

Increased  %  hole  area  "Ac"  in  the  cage  worsens  across-the-wafer 
uniformity  (Fig. 5. 3)  because  of  the  rise  in  the  amount  of  SiH:  flowing 
in  from  the  annulus.  Batch  variance  (Fig.  5.4)  is  not  influenced 
strongly  by  "Ac"  for  lower  levels  but  at  20%  hole  area  the 
consumption  of  the  annular  SiH2  is  so  rapid  that  the  increased 
depletion  effects  lead  to  a  rise  in  batch  variance.  Therefore,  SNt2 
falls  rapidly.  Again,  the  increased  consumption  of  the  annular  SiH2 
leads  to  slightly  increased  average  growth  rates  (Fig. 5. 5). 
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Fig.  5.4:  Marginal  Plots  for  SNt2 

The  film  thickness  uniformity  degrades  and  the  average  growth  rate 
rises  with  increasing  deposition  temperature  "Td".  Increasing  growth 
rate  with  temperature  is  explained  by  the  rising  rate  of  production  of 
SiH:  (the  dominant  surface  reactant)  in  the  homogeneous  phase 
(Fig. 5. 5).  However,  film  nonuniformity  rises  for  temperature  upto 
600^C  and  then  changes  very  little  (Fig. 5. 3  and  Fig.  5.4).  Probable 
cause  for  this  behavior  is  the  relative  contribution  of  various  gas 
phase  species  involved  in  the  polysilicon  film  deposition.  Upto  600°C 
the  increased  contribution  of  SiH:  from  the  larger  volume  of  the 
annular  region  accounts  for  the  nonuniformity.  For  higher 
temperatures  other  highly  reactive  gas  phase  minority  species  like 
disilane  (Si:H6).  that  add  to  the  nonuniformity,  become  extremely 
unstable  and.  therefore,  their  contribution  to  the  heterogeneous 
chemistry  (i.e  surface  reactions)  is  reduced.  Also,  at  nigher 
temperature  silane  (SiH4;,  the  majority  species  in  the  gas  phase, 
becomes  more  active  in  attacking  the  surface  sites  successfully  and 
this  makes  the  growth  rate  kinetics  tend  towards  kinetic  control  and 
higher  uniformity  rather  than  transport  limited  control.  Therefore, 
there  is  no  further  rise  in  the  nonuniformity.  The  existence  of  these 
various  parallel  reaction  mechanisms  has  been  well  documented  [4, 
5]  but  the  experimental  verification  of  the  above  reasoning  is 
required. 
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Fig.  5.5:  Marginal  Plots  for  SNg 

The  influence  of  rising  pressure  on  batch  variance  (Fig.  5.4)  is  not 
verv  significant  but  the  film  nonuniformitv  across  the  wafer  rises 
appreciably  with  rising  pressure  (Fig. 5. 3).  Rising  pressure  means  a 
higher  SiH;  concentration  and,  therefore,  a  higher  growth  rate 
(Fig. 5. 5 1.  .Although  rising  pressure  means  a  higher  phosphine  (PH?) 
concentration  as  well,  it  can  be  shown  by  order  of  magnitude 
analysis  that  the  change  in  the  PH3  concentration  with  pressure  is 
louer  than  that  of  SiH2.  In  any  case,  the  rise  in  PH3  concentration 
does  have  a  dampening  effect  on  the  rise  in  growth  rate  and. 
therefore,  the  rise  shown  in  Fig. 5. 5  is  not  very  significant. 

The  overall  batch  variance  improves  for  the  middle  level  of  inier- 
injector  distance  "lx"  (Fig.  5.4)  but  the  across-the-wafer  variance  and 
growth  rate  are  worst  for  this  level  (Fig. 5. 3  and  5.5).  The  influence 
of  "Ix"  on  the  film  thickness  properties  is  not  well  understood  except 
for  a  possible  linkage  between  gas  phase  SiH3PH2  complex  and  the 
point  where  the  gases  from  injector  #2  are  released  into  the  tube 
u'ith  respect  to  the  position  of  the  heated  load  zone.  Th  ^'  is  briefly 
discussed  in  section  3.2  of  chapter  3. 


3  7 


Thickness  variance  invariably  falls  (Fig5.3  and  Fig.  5.4)  and  growth 
rate  rises  (Fig. 5. 5)  with  rising  silane  flow  rates.  Uniformity  across 
the  wafer  improves  with  rising  silane  flow  rates  possibly  because 
part  of  annular  SiH2  is  blown  away  at  higher  silane  flows. 

As  e.xpected,  the  thickness  variance  across  the  wafer  rises  with  an 
increase  in  the  phosphine  flow  rates  (Fig. 5. 3  and  Fig.  5.4).  This  is 
because  of  the  change  in  the  reaction  mechanism  from  being  kinetic 
limited  to  transport  limited  with  rising  dopant  gas  concentration. 
Increased  PH3  means  more  poisoning  of  the  substrate  and,  therefore, 
a  substantially  lower  growth  rate(Fig.5.5).  The  influence  of  PH3  on 
batch  variance  is  different  for  "Qpi"  and  "Qp2"  probably  because  of 
the  effect  of  "Ix". 


5.2  Resistivity  Data  Analysis: 

Marginal  plots  in  Fig. 5. 6  show  the  main  effects  of  the  process 
parameters  on  the  average  resistivity  of  the  in-situ  phosphorus 
doped  polysilicon  film.  The  corresponding  ANOVA  analysis  for 
smaller-the-better  resistivity  SN  ratio  (SNp)  is  given  in  Table  A. 5  in 
Appendi.x  A.  Phosphine  flow'  rates,  interwafer  distance,  inter- 


Fig.  5.6:  Marginal  Plots  for  SNp 
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injector  distance,  and  cage  hole  area  have  major  effect  on  resistivity. 
Note  that  the  effort  here  was  to  reduce  the  average  resistivity  across 
the  wafer  in  a  load. 

Increased  "w"  means  a  higher  growth  rate  of  the  film  and,  therefore, 
larger  grain  size.  Also,  resistivity  is  seen  to  fall  with  rising  thickness 
[9],  Therefore,  SNp  rises  with  increase  in  "w". 

Increasing  film  growth  rate  with  rising  "Ac"  leads  to  a  higher 

resistivity  (i.e.  lower  ^Np)  but  the  effect  reverses  for  very  high  "Ac". 
It  is  also  observed  that  the  variance  of  resistivity  across  a  wafer 
increases  dramatically  for  high  levels  of  "Ac".  This  indicates  that 
there  is  an  increased  contribution  of  PH3  from  the  annular  region  to 
the  interwafer  region  at  very  high  hole  area.  Therefore,  for  large 

holes,  the  probable  cause  for  lower  resistivity  is  that  the  increased 

influx  of  annular  PH3  content  overshadows  the  rise  in  thickness  due 

to  increasing  influx  of  annular  SiH2  and  this  raises  the  dopant 
densit}'  in  the  film. 

Deposition  temperature  rise  leads  to  higher  resistivity  because  of  the 
enhanced  desorption  of  phosphorus  from  the  silicon  surface.  At 
higher  deposition  temperature  the  grain  size  of  polysilicon  is  small 
and  this  too  leads  to  higher  resistivity. 

Rising  gas  pressure  means  an  increased  partial  pressure  of  PH3.  This 
means  a  rise  in  the  phosphorus  content  of  the  film  and,  therefore,  a 
fall  in  resistivity. 

There  is  dramatic  rise  in  resistivity  with  increase  in  the  inter¬ 
injector  distance.  This  behavior  can  again  be  explained  by  the  role  of 
gas  phase  SiH3PH2  complex.  As  discussed  in  Chapter  3,  at  highest 
level  of  Ix  more  of  phosphine  from  injector  #2  is  trapped  in  a 
complex.  Therefore,  lesser  amount  of  phosphine  is  available  for  the 
wafers.  This  leads  to  lower  dopant  density  in  the  film  and,  hence, 
higher  resistivity. 
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It  has  already  been  shown  that  a  rise  in  silane  flow  leads  to  thicker 
film.  Also,  increased  silane  flow  rate  blows  aw-ay  part  of  the  annular 
PH  3  (this  is  indicated  by  reduced  variance  in  resistivity  that  is 
observed  across  the  wafer  with  rising  silane  flow').  These  two  facts 
explain  the  observation  that  there  is  a  rise  in  average  resistivity  with 
increased  silane  flow. 

Increased  PH3  flows  increase  the  concentration  of  phosphine  in  the 
system.  This  leads  to  greater  amount  of  phosphorus  incorporation 
into  the  polysilicon  film  and,  therefore,  lower  resistivity. 


5.3  Verification  Experiments: 

In  order  to  study  the  effectiveness  of  the  orthogonal  design  of 
experiments,  a  set  of  confirmation  runs  were  made  based  on  the 
optimum  levels  for  the  process  parameters  predicted  by  the 
marginal  plots.  These  optimum  levels  correspond  to  the  highest  SN 
ratio  in  the  relevant  marginal  plots  in  Fig. 5. 3,  5.4,  5.5  &  5.6.  The  data 
for  these  runs  is  given  in  Table  A. 6  in  Appendix  A.  Two 
modifications  were  made  to  the  normal  procedure  of  determining  the 
optimum  level  of  a  process  parameter  as  the  one  with  highest  SN 
ratio.  Firstly,  in  order  to  avoid  cutting  the  load  size,  only  lower  level 
of  interwafer  distance  was  used  for  the  confirmation  runs.  Secondly, 
instead  of  using  the  optimum  levels  from  the  marginal  graphs, 
second  order,  16  terms  linear  regression  equations  for  various  SN 
ratios  were  evolved  by  using  the  18  data  points  for  each  of  the  SN 
ratios.  For  instance,  one  of  the  regression  equations  was^: 


'  The  units  various  parameters  in  equation  5.1  are: 

[wj  =[-0-1  i.c.  just  level  #:  [Ad  =  [TdJ  =  pC);  [P]  =  [mtorr];  [lx]  =  [inches]; 
[Qsl  =  [seem];  [Qpd  =  [Qp:I  =  [seem] 
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+ 


SNp=  2.367E+02  +  (2.919E+00)(w)  -  (8.387E-01)(AJ  +  (4.197E-02)(aJ) 
(5.570E-01)(Tp)  +  (4.360E-04)(T‘)  -  (6.982E-01)(P)  +  (1.047E-03)(p2) 
(1.175E-K)0XIJ  -  (3.654E-02)(l;)  +  (1.279E-01)(Q,)  - 
{2A3OE-O40  +(4.268E+00)(Qpj)  -  (5.556E-01)lQ2pJ  + 
(5.447E+00)(Qp2)  -  A53E+OOi(^J 


From  such  equations  the  extrema  regression  predictions  for  the 
maxima  of  the  relevant  SN  ratio  were  determined  as  the  optimum. 


SN 

Type 

Process 

Parameters  for 
confirmation 

Runs 

SN  ratio  Values 

Response  Back-Calculated 
from  Corresponding  SN 

predicted 

observed 

Response 

Type 

Value 

Confirm¬ 
ation  Run 

S^, 

w  =  .75;  =  .125; 

T,  =  853  .;  P  =  325  .; 
1,  =  65.;  Os  =  350.2; 
0p,=  .58,0.,=  .26 

49.4 

46.2 

/4cyv>v  /W 
HK/Ziv 
!•  WJjV 

flhnonfrwniffv 

0.45% 

(run  *1) 

0.49% 

SN,, 

V  =  .75;  =  .328; 

=  853.,  P  =  375.; 
I,=  53.5;  Os  »  405.9; 
Opr  2.87;  0*,=  .79 

31.1 

32.2 

MMHjnifwnirfy 
rb/-  j 

(M/vsy  ii'jAv-  1* 
uOii')7  /IV 

4.7  % 

(run  *2) 

2.5  % 

SN3 

v=  .75;  r„  =  .783, 

=  893.,P  =  375.; 
I,  =  65.;  Os  *350.2, 

Ctr  -58.  V  -26 

42.2 

40.7 

iJfVirfli  ^J/V 
Av  J  SjfcA 

61.2 

A^/min 

(run  *5) 

98.0 

A*Vmin 

SNp 

v=  .75,r,  =  .125; 

T,  =  853.;  P  =  375.; 
I,=  42.2;  Os*  264.2, 
Opr  3.07;  0^,,*  1  -52 

-22.3 

-21.9 

rbf'j  Bijfcfi 

12.45e-4 

ohm-cm 

Note:  Units  for  Process  Peremeters  are  Iv]=[r  ]=llj  =Icm.];  [Tjj]=l®K]; 
[Pl=[nntorr];  (QsJ  =[Qpil*[C^2^ 


Table  5.3:  Results  of  Confirmation  Runs  vs  L18  Orthogonal 
Arrau  Experimentation 


For  example,  in  equation  5.1  the  derivative  of  SN’  with  respect  to  "Ix" 
is  equated  equal  to  zero  to  get  a  predicted  optimum  of  Ix  =  16.6 
inches.  In  case  of  process  parameters  with  negative  slopes  the 
extrema  regression  prediction  is  for  the  minima  of  the  concerned  SN 
ratio.  Therefore,  the  marginal  graph  predictions  for  the  optimum  (i.e. 
maximum  SN)  were  used  for  parameters  with  negative  slopes.  For 
example,  the  extrema  for  "Ac"  corresponds  to  a  minima  in  equation 
5.1.  The  process  parameter  settings  and  the  predicted  SN  ratios  for 
the  confirmation  experiments  are  shown  in  Table  5.3. 

As  indicated  before,  a  lack  of  data  for  determining  the  cost 
coefficients  in  equation  2.11  prevented  us  from  performing  a 
confirmation  experiment  for  the  overall  optimum  for  the  tradeoff 
responses  of  film  thickness  uniformity,  growth  rate,  and  resistivity. 

Results  of  these  confirmation  experiments  are  tabulated  in  Table  5.3 
which  shows  the  observed  and  predicted  SN  ratios  for  these  runs. 
The  predicted  SN  ratios  seem  to  be  right  on  target  except  for  SNji. 
This  indicates  that  the  additive  model  underlying  the  Taguchi 
method  of  process  optimization  is  not  sufficient  for  film  thickness 
uniformity  across  the  wafer.  There  are  some  important  interactions 
that  need  to  be  considered  for  determining  the  optimal  set  of  process 
conditions.  Development  of  a  semi-empirical  model  serves  to  this 
end. 

The  confirmation  run  result  for  SNp  is  not  as  good  as  the  best  run  in 
LI 8  orthogonal  array  (run  #17)  because  for  the  actual  optimal  run  it 
was  required  that  highest  levels  of  interwafer  distance  and  cage  hole 
area  be  used  for  which  the  predicted  optimal  SNp  y^as  -8.9  i.e.  a 
resistivity  of  .28  mohm-cm.  But  we  used  lowest  level  of  "w"  for  the 
reason  stated  above.  Also,  a  cage  hole  area  of  0.5%  was  used  so  that 
the  average  resistivity  is  not  only  reasonably  low  (which  1.245 
mohm-cm  is)  but  also  uniform  across  the  wafer.  It  is  important  to 
point  out  that  the  average  nonuniformity  for  sheet  resistance  in  run 
#17  was  11%  wheras  for  the  SNp  confirmation  run  it  was  only  4%. 


Regression 

Coefficients 

Primitive  Variables  Model 

Dimensional  Analysis  Model 

_ Estimated  Values _ 

_ Estimated  Values 

Wafer#2 

Wafer#7 

Wafer#l  l 

Wafer«2 

Wafer#7 

Wafer#1 1 

0) 

1.150e+2 

1.267e+2 

1.71 1e+2 

8.072e+1 

7.084e+1 

1 .228e+2 

Cl 

2.0316-1 

7.004e-2 

-1.994e+0 

-1.695e-t-( 

1  -2.266e+C 

-4.31 8e+0 

C2 

■4,278e+0 

-5.466e+0 

-4.162e+0 

-7.427e+0 

-9.3706+0 

-8.747e+C 

C3 

-2.302e-2 

-2.831  e-3 

-2.978e-2 

-6.572e-2 

-4.0648-2 

-9.004e-2 

C4 

-3.945e-2 

-5.695e-2 

-7.789e-2 

-8.089e-1 

-1.958e+0 

-2.826e+C 

C5 

-3.639e+0 

-5.320e+0 

-6.246e-*.0 

-1.616e-1 

-2.882e-1 

-1.698e-1 

Cg 

1.498e-2 

5.297e-2 

-2.247e-2 

1.153e-1 

2.524e-1 

1.177e-1 

C7 

-5.103e-2  • 

1 .388e+0 

-1.043e+0 

1.214e-2 

1 .733e-2 

2.086e-2 

Cg 

4,148e-2 

-2.645e+0 

-2.190e+0 

-2.6608-: 

-4.373e-5 

-4.280e-5 

Cg 

2.149d-2 

5.33 9e-2 

7.254e-2 

2.092e-2 

5.214e-2 

7.362e-2 

Cl  0 

-2.367e-4 

-2.999e-4 

-2.334e-4 

-1.282e-4 

-2.808e-4 

-1.319e-4 

C11 

9.705e-2 

1.053e-2 

1.182e-2 

1.232e-3 

1 .564e-3 

2.067e-3 

O2 

-4.0548-2 

2.200e-2 

-2.290e-2 

-5.160e-4 

-6.908e-4 

-8.243e--i 

Table  5.4:  Calibrated  Repression  Coefficients  for  Primitive 
Variable  Regression  and  Dimensional  Analysis  Models 


5.4  Model  Calibration  and  Analysis 

The  dimensional  analysis  model  evolved  in  chapter  3  for  the  SN  ratio 
(defined  in  equation  3.1)  was  regressed  to  the  thickness  mean  and 

variance  data  .  Instead  of  developing  a  single  generalized  SN  ratio 
model  for  a  batch,  the  data  was  utilized  to  develop  a  model  for 

individual  wafers  located  at  different  axial  positions  in  the  load.  This 
approach  did  not  require  any  extra  effort  in  terms  of 

experimentation  or  data  acquisition.  In  fact,  by  doing  so,  the 

accuracy  of  the  calibrated  models  was  not  degraded  by  the  lack  of  fit 
for  axial  data  that  would  otherwise  be  introduced  if  the  axial 


w 

A 

F 

Primitive 

Variables  Model 

Dimer)sional 
Analysis  Model 

E 

R 

it 

F-Test 

u. - 

O 

> 

< 

F-Test 

Av.  %  Error 

2 

4.1 

5.8 

27.9 

2.4 

3 

3.8 

6.4 

48.0 

2.0 

4 

4.2 

6.2 

35.0 

2.1 

5 

4.1 

7.1 

41.7 

2.6 

6 

5.3 

6.6 

49.6 

2.1 

7 

5.7 

6.5 

47.3 

2.0 

8 

5.5 

6.1 

19.7 

3.2 

9 

5.8 

6.5 

45.3 

2.3 

1  0 

6.8 

5.7 

37.6 

2.5 

1  1 

7.8 

5.7 

46.8 

2.5 

Table  5.5:  F-Test  &  Average  %Error  Comparison 
of  Primitive  Model  vs.  Dimensional  Model 


dependence  was  considered.  Therefore,  coefficients  Co,  Cj .  and 

C 12  in  equation  3.12  were  determined  for  SN  models  for  all  test 
wafers  in  each  run  by  multiple  linear  regression.  The  calibrated 
values  of  these  coefficients  for  test  wafers  #2,  #7,  and  #11  (see 
Fig.4.1)  are  given  in  Table  5.4. 

For  the  same  sets  of  data  points,  the  SN  model  in  equation  3.13  was 
calibrated.  The  calibrated  values  for  the  regression  coefficients  are 
given  in  Table  5.4. 

Table  5.5  deals  with  the  comparative  analysis  of  the  models  in 
equation  3.12  and  3.13.  It  is  very  clear  that  the  semi-empirical 
model  based  on  dimensional  analysis  is  much  superior  than  the 
corresponding  primitive  variables  model.  The  former  has  a  much 
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lower  average  %error^  and  a  F-ratio^  that  is  six  to  seven  times  better 
than  the  latter. 

The  accuracy  and  better  predictive  power  of  the  semi-empirical 
models  is  further  demonstrated  by  predicting  SN  ratios  for  the 
individual  wafers  for  certain  sets  of  process  parameters.  These 


TESTING  CONDITIONS 

Wafer 

» 

§N 

obser 

-ved 

SN 

predicted. 

Primitive 

Run  # 

Process  Parameters 

1. 

w  =  .75  cm:  r^,  =  .125  cm; 

To  =  853.  K:  P  =  375.  mtorr; 

1.  =  49.8  cm;  Cfe  =  288.3  seem; 
Q-  =  3.14sccm:  Qjz  =  1.07sccm 

2 

28.4 

29.2  +/-  7.6 

a 

31.2 

28.4  +/-  8.2 

30.6  +/-  2.0 

1 1 

32.6 

30.8  +/•  8-6 

34.1  +/-  2.6 

2. 

w  =  .75  cm;  r-  =  .125  cm; 

To  =  853.  K;  P  =  325.  mtorr; 
li  =  65.1  cm;  Cfe  =  350.2  seem; 
Q,-  =  .58  seem:  Q)2=  26  seem 

2 

44.5 

47.4+/-11.4 

47.5  +/-  4.1 

n 

49.0 

52.6+/-12.2 

54.5  +/-  4.0 

1 1 

.... 

65.6+/-13.0 

66.8  +/-  5.0 

3. 

w  =  .75  cm;  r^  =  .125  cm; 

To  =  853.  K;  P  =  375.  mtorr; 

1,  =  42.2  cm;  =  264.2  seem; 

Q,  =  3.07sccm;  Q,2  =  1.52sccm 

2 

29.0 

31.4+/-13.4 

30.6  +/-  2.6 

|[Q|| 

33.6 

29.4+/-14.4 

30.6  +/-  2.4 

12 

30.2 

26.0  +/■  8.6 

29.4  +/-  4.1 

Table  5.6:  SN  Ratios  Observed  vs  Predicted  bv  Primitive 
Regression,  and  Dimensional  Analysis  Models 


results  for  selected  wafers  are  shown  in  Table  5.6  (results  do  hold  for 
all  the  wafers  in  a  batch)3.  The  dimensional  analysis  model  predicts 
as  good  as  the  primitive  variables  model  but  the  predictions  of 
dimensional  analysis  model  are  much  more  reliable  because  of  their 
narrow  confidence  intervals.  This  high  reliability  of  the  dimensional 
analysis  model  stems  from  the  incorporation  of  the  process  physics 
into  the  model. 


'  The  average  %error  is  defined  as  the  average  of  the  %  absolute  deviation  of 
the  predicted  SN  value  from  the  observed  value. 

2  The  F-ratio  is  the  ratio  of  the  mean  square  regression  of  calibartion 
coefficients  to  the  mean  square  error  in  the  regression. 

3  Note  that  the  SN  ratio  in  Table  5.6  is  for  individual  wafers  in  a  batch  and  not 
the  SN  ratio  as  used  in  Taguchi's  method  before. 
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The  semi-empirical  models  for  individual  wafers  were  used  to 
predict  the  batch's  (and  not  individual  wafers’s)  SN  ratios  SNxi  and 
SN't:  for  the  confirmation  runs  and  the  validity  of  these  predictions 
was  compared  to  that  of  Taguchi’s  averaging  method.  The  SNji  from 
dimensional  analysis  models,  called  (SNyOd,  model*  determined  by 
back-calculating  ratio  for  every  i'*’  wafer  from  its 

dimensional  analysis  model.  These  ratios  for  all  wafers  in  a  batch 
were  added  and  averaged  and  then  expressed  in  a  decibel  scale  (as 
in  equation  2.2  for  SNxi)  to  give  the  value  for  (SNxi)d,  model- 
Similarly,  the  SNxi  from  primitive  variables  models,  called 
( SNxi  )p, model-  ^as  determined  from  the  primitive  variables  models 
for  individual  wafers.  The  SNx2's  were  determined  by  regressing 
equations  3.12  and  3.13  respectively  to  the  data  values  of  SNx2  for 
18  runs.  This  gave  semi-empirical  model  equations  for  (SNx2)d,  model 
(from  dimensional  analysis  equation  3.12)  and  (SNx2)p,  model  (from 
primitive  variables  equation  3.13).  These  equations  for  (SNx2)d,  model 
and  (SNx:)p,  model  '^'sre  used  for  predicting  the  SNx2  for  the 
corresponding  confirmation  run  conditions.  All  these  predictions  in 
comparison  to  those  of  Taguchi's  averaging  method  are  shown  in 


Table  5.7. 


SN 

Type 

Process 

Parameters  tor 
confirmation 

Runs 

SN  Ratio  Values  for  Confirmation  Runs 

observed 

predicted, 

Taguchi's 

predicted. 

Dimensional 

predicted, 

Primitive 

SN,, 

w  =  .75;  r,  =  .125; 

Tj  =  853.;  P  =  325.; 
I,  =  65.;  Q  =  350.2; 
a-  =  .58;  Q2=  .26 

46.2 

49.4 
+  /- 

3.0 

47.2 

+/- 

3.7 

42.4 

+/- 

10.2 

1 

w  =  .75;  r„  =  .328; 

Ta  =  853.;  P  =  375.; 

I,  =  53.5;  Q  =  405.9; 
Q,  =  2.87;  0,2=  .79 

32.2 

B 

30.3 

+  /  ■ 

3.3 

B3B 

Note:  Units  for  Process  Parameters  are:  (w]=[rh)=[lx]=[cm.]:  [Td]=[°K]: 


[P]=[mtorr]:  [Qsl  =  [Qpi  ]  =  [Qp2]  =  [sccm) 

Table  5.7:  Model  Predictions  vs.  Taauchi's  Method 


Once  again  the  predictions  by  dimensional  analysis  models  are  more 
reliable  than  those  by  primitive  variable  models  because  the  former 
has  a  narrower  confidence  interval.  Also,  these  model  predictions 
are  better  than  those  by  Taguchi’s  averaging  method  in  terms  of 
confidence  interval  for  SNt2  and  mean  value  for  SNji.  This 
improvement  is  accounted  for  by  the  fact  that  in  evolving  these 
semi-empirical  models  we  have  incorporated  important  process 
variable  interactions  that  Taguchi’s  method  tends  to  neglect. 


6.  CONCLUSIONS 


The  orthogonal  array  design  of  experiments  for  process  and 
equipment  optimization  was  combined  with  a  semi-empirical 
modelling  approach  based  on  dimensional  analysis.  Typical 
constraints  of  a  purely  empirical  approach  were  experienced  while 
applying  Taguchi's  method.  These  constraints  included  the  lack  of 
access  to  the  desired  experimental  environment,  a  general  drawback 
in  empirical  approach,  and  the  problem  of  process  parameter 
interactions,  a  particular  limitation  to  Taguchi's  orthogonal  array 
method.  Taguchi's  method  worked  well  for  the  optimization  of  In- 
Situ  Phosphorus  Doped  Polysilicon  Process  via  LPCVD.  However,  its 
dependence  on  the  assumption  that  the  process  under  consideration 
is  free  of  parameter  interactions  proved  to  be  an  inhibiting  limitation 
for  at  least  one  of  the  responses  studied  i.e  nominal-is-the-best 
signal-to-noise  ratio  SNji  for  film  thickness  uniformity  across  the 
wafer.  As  a  solution  to  such  problems  of  purely  experimental 
approach  it  was  recommended  to  evolve  a  semi-empirical  model  for 
the  process  under  consideration.  Dimensional  analysis  was  used  as  a 
way  of  formulating  such  a  model  for  the  film  thickness  uniformity  in 
the  process  under  consideration. 

Dimensional  analysis  was  used  to  transform  primitive  variables  that 
described  the  process  under  consideration  into  dimensionless  groups. 
It  was  theorized  that  these  dimensionless  groups  would  enhance 
empirical  modeling  by  infusing  some  physical  or  causal  reality  of  the 
process  into  the  model.  Enhancement,  as  expected,  occured  on  two 
fronts.  Dimensional  analysis  reduced  the  number  of  primitive 
experimental  variables  from  eight  to  seven  and  primitive  process 
variables  (of  which  only  8  were  varied  in  the  experiments)  from 
twelve  to  seven.  Secondly,  dimensional  analysis,  as  theorized, 
yielded  dimensionless  parameters  that  were  physically  more 
significant  than  the  primitive  variables.  The  average  %  error  in 
dimensional  analysis  model  was  about  a  third  of  that  in  primitive 


48 


variables  model.  Also,  from  statistical  point  of  view,  the  F-test  for 
dimensional  analysis  model  was,  on  the  average,  six  to  seven  times 
better  than  that  for  the  corresponding  primitive  variables  model. 
The  dimensional  analysis  models  had  a  more  accurate  and  reliable 
predictive  power  as  measured  by  the  fact  that  the  confidence 
interval  of  its  predictions  for  a  certain  set  of  process  conditions  was 
more  than  three  times  narrower  than  that  for  primitive  variables 
model  predictions.  The  semi-empirical  model  overcame  the  problem 
encountered  by  Taguchi's  method  with  regards  to  the  importance  of 
the  process  parameter  interactions  that  the  later  tends  to  neglect. 

For  future  work  it  is  recommended  that  the  semi-empirical  model  be 
derived  for  a  response  corresponding  to  the  SNoveraii  defined  before. 
Dimensional  analysis,  as  demonstrated  in  this  work,  should  be  used 
to  evolve  a  model  for  such  overall  system  response.  For  this  to 
happen  it  is  essential  that  the  desired  cost  coefficients  be  determined 
by  studying  the  economics  of  semi-conductor  manufacturing.  Such  a 
model  should  be  used  to  design  a  D-optimal  set  of  experiments  for 
optimizing  and.  therefore,  designing  process  and  equipment  for  in- 
situ  doped  polysilicon  LPCVD  process.  It  is  also  recommended  that 
the  strength  of  the  semi-empirical  approach  shown  in  this  work 

should  be  tested  against  a  purely  mechanistic  model  for  the  same 

process.  Such  a  comparison  is  expected  to  reinforce  the  validity  of 
the  point  made  in  this  work  that  semi-empirical  approach  is  the 
optimal  way  of  designing  optimal  processes  and  equipment  from 

modelling  point  of  view.  The  potential  of  such  an  approach  to 

process  modelling  in  Manufacturing  is  enormous. 
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APPENDIX  A:  Data  Tables  and  Plots 


The  data  acquired  for  the  film  thickness  and  resistivity  of  all  the 
wafers  in  each  of  the  18  runs  is  outlined  below  in  Tables  A.l(l)  to 
A.  1(18).  The  .ANOV.A  analysis  for  these  runs  is  shown  in  Tables  A. 2, 
.A. 3,  A. 4,  and  A. 5.  The  data  for  confirmatin  runs  is  shown  in  Table 
A. 6. 

In  the  data  Tables  A. 1(1)  to  A. 1(18)  and  Table  A. 6,  ”t"  is  the  film 
thickness.  "Std.Dev._t"  is  the  ^  standard  deviation  in  film  thickness, 
"rho"  refers  to  the  resisti\ity  of  the  film,  and  "Std.Dev._rho"  is  the 
estimated'  standard  deviation  of  resistivity.  Next  to  the  data 
tables,  the  axial  profiles  of  average  film  thickness  and  resistivity  are 
shown  for  each  run.  The  missing  data  in  table  A. 6  is  due  to  the  fact 
that  the  availability  of  insufficient  number  of  test  wafers  forced  us  to 
replace  the  test  wafers  in  some  of  the  slots  in  the  load  by  the  dummy 
wafers.  Also,  some  of  the  test  wafers  turned  out  to  be  defective  and. 
therefore,  accounted  for  the  missing  data  points. 


'  Standard  deviation  in  resistivity  is  estimated  by  multiplsmg  standard 
deviations  of  Him  thickness  and  sheet  resistance. 


Run  #1 


t  (  Args::rm 

3340  COO 
3462.500 
3533. 1  00 
3691.600 
3789  600 
3914  SCO 
4075  800 
4286.300 
4552  ’  00 
4736.800 
4669  ’  00 
3856. 103 


data. rune  1 

5000 

4800 

Std.  Cev  _t 

rtio  (iff JT—Ctii) 

Sld.Dev._rho 

fl 

_46C0 

O 

gl400 

1 .07"i 

8  852 

0.021  r. 

c 

A 

5;200 

0.54% 

0  33% 

0  32”. 

8.248 

8  681 

9  281 

0.004% 

0.004".'. 

0.001 

U 

X 

E 

^’.000 
5  1800 
~I6C0 

0,23”. 

0.21”'. 

10  077 

1  0.992 

0  002”'. 
0.002”'. 

(Z 

1400 

0.23”'. 

12  293 

0  0017'. 

3200 

0,36”'. 

13.986 

0  003“'. 

3000 

16  237 
18.417 
16  790 
13.940 


0  0C2-. 
Q  001% 
0  00 1-. 
0,008% 


2  4  8  8 

Wifar  Number 


8  10  12 


data  run»2 


Run  M2 


1  (  Angsiorrr  ) 

2  25'5, 00000 

3  2550  00000 

4  2553  80000 

5  2523  30000 

3  2433  50000 
7  2487  '  0003 

3  2550  50000 
9  2543  80000 

*0  2'’70’0000 
•  '  29-  '  80000 

12  2'12  80003 


Sid 

D9v,_’ 

rr-o  (10'.A_cm) 

Sid.Dav._rno 

3600 

3400 

4 

87”. 

6  885 

0  2817'. 

3200 

2 

8  5”. 

7. 73' 

0  064”'. 

n 

in 

'“1000 

2 

2 

5  0”. 

3’  ”, 
S'", 

7  32' 
7.905 

7  964 

0  075*; 
0.082% 
0.068"; 

c 

u 

X 

K 

1 1300 

|;600 

tr, 

2 

25"', 

8  043 

0  057-, 

:  '400 

.f 

2 

3  2"'. 

8,2'2 

0  058"; 

E 

:2oo 

2 

20", 

8  42T 

0  058”; 

ul 

!000 

1 

55". 

8  539 

0.027”; 

1800 

1 

53", 

8  '5' 

0  02'”. 

1600 

I 

36", 

8  '90 

0  019”; 

2 

33"', 

7  88" 

0.060”; 

2  4  6  81012 

Water  Number 


Ca:a.runa3 


1(A) 

Sid  Car  _! 

mo  ;i07i__c-ri 

Std  Dav  _rrio 

3800 

1 

32"*^ 

7  7  2”; 

8  854 

0  295”'. 

3600 

2 

2353  00000 

5  5  3". 

8  3-1 

0  151”; 

„  3400 

3 

2'33  50000 

4  35”, 

8  67  1 

0  144”', 

S  -1200 

4 

2552  80000 

5  71"'. 

8  582 

0  239”; 

i  SiCCO 

5 

2455  '0000 

5  51”  ', 

8  543 

0  212”; 

u  2 

2  ^  ’300 

6 

2373  10000 

5  20", 

8  S'6 

0  212”; 

*”  c’l600 

7 

2235  OCCOO 

4  68". 

8  68" 

0  178”; 

E  —’400 

8 

2258  00000 

3  5  5”; 

8  8-5 

0  098”; 

\L  !200 

9 

2325  30000 

3  95"; 

8  986 

0  1257; 

1  0 

2458  30000 

4  91”; 

9,192 

0.191”; 

fcUUlJ 

1  1 

2473  10000 

3  90”', 

9,128 

0  113”; 

1 630 

1  2 

25  I  8.30000 

6,63”.; 

8.132 

0.331% 

Run  #3 


2  4  6  8  10  12 

Wafer  Number 


Table  A.KD- 


Data  from  Runs  1.  2.  and  3 


Realsllvlty  (1o  ''_n.  cm)  ReslatIvMy  (10  _cm)  Realallvlly  (I0  ‘'jz._cm) 


(>  <I) 


Run  #4 


1  (Angs'.orm) 

Sid. Dev. _i 

rtio  (10'2i_cm) 

1 

1038.10000 

7.20". 

9  148 

2 

1  145.80000 

5.04“'. 

8.729 

3 

1  152.50000 

5.49“. 

8  779 

4 

1  1  07.50000 

5.38“'. 

8  902 

5 

1  054.00000 

5.92“'. 

9  207 

6 

1  009.00000 

6  05“'. 

9.398 

7 

970,1  6000 

6.66". 

9  529 

8 

938  85000 

7.09“'. 

9.498 

9 

937.56000 

5.3  9". 

9.644 

0 

1 021.90000 

4.95"'. 

9.536 

1  009.00000 

4  53“'. 

9,510 

2 

866.23000 

5.91“'. 

9.88' 

data.run«4 

Std.D*v._fho 

O.SOC% 
0,254% 
0.359% 
0.260% 
0  416% 
0.475% 
0.545% 
0  568% 
0.335?'. 
0  284?'. 
0.242% 
0.498?'. 


11.3 
10  9  i 
10.5  ", 
10.1 
9.7 

9.3  C 
8.9  > 


7.3 


Wafer  Number 


daia. runes 


(ArgsiCfT-.i 

Sic. Dev. _i 

rtio  (13”r_cm) 

Sid.Dav._rho 

^3CS 

2  1  0“'. 

10.842 

0.047“'. 

65'9  SOOOC 

1 .'  9“. 

14  322 

0  030“'. 

6893  70000 

C  S3’. 

18,222 

0  018“'. 

706'  80000 

C  8*"'. 

2'  754 

0.018?; 

7-  65.20000 

0  3  0“. 

25  436 

0  013“'. 

7253  80000 

0  Sc“'. 

23  246 

c  011“; 

7354  000:0 

0  82“'. 

30  9-2 

0  012“. 

*552  30000 

0  82“'. 

35  750 

0  013“'. 

8033  SOOCO 

0  38“. 

45  520 

0.007% 

8483  40000 

0  25“. 

56  850 

0  005“'. 

34*3 

0  16“. 

64  322 

0  003“'. 

7137  30000 

0,51“. 

50  989 

0  010". 

Run  #5 


70 


60 


50 


40  Z, 


20  •= 


10 


Wafer  Number 


da'.a.runee 


1  (Argsiorm  )  Sld-Dev_l  rtio  (10'J;_crn)  Sid  Dev.  rho 


1  36' 2 

2  38-3 

3  3670 

4  3499 

5  3376 

6  3281 

7  3252 

8  3324 

9  34-8 
10  3551 
1  1  3593 
12  3254 


60000 
30000 
30000 
50000 
80000 
80000 
30000 
1  0000 
60000 
60000 
60000 
60000 


57' 
89“ 
1  5“ 


6 

5 

6 

5  1  7“ 
5.77“ 
4  96“ 
4  2'“ 
4  25“ 
3  99“ 
3  35“ 
2.75" 
3.70“ 


1  1  250 
10  989 

1  0  49' 

10  191 
10  012 
9  885 
9  790 
9  909 
10  297 
1  0  67' 
11129 
1  0  698 


0  22'“'. 
0  184“', 
0  238“'. 
0.168“i 
0.21  6?'. 
0  171?'. 
0  139". 
0  154?'. 
0  109*'. 
0  104". 
0  060". 
0,1 18?'. 


Run  #6 


12.5 


11.7 


10 : 
10.1 
9,7 
9.3 
8  9 
8.5 


Table  A.l(4)'(6): 


Data  from  Runs  4.  5.  and  6 


data.run#? 


Run  U7 


t  (Angsiorm)  Stfi.Dev.j  mo  (10'^_cm)  Sld.D«v._rho 


1  2334.60000 

2  23S4  30000 

3  2916.30000 

4  2369.80000 

5  2679.00000 

6  2259  3QC00 

7  1687.60000 

8  1623.00000 

9  1605.40000 


0.234% 

0.269% 

0.174% 

0.186% 

0.188% 

0.741% 

0.789% 

0.674% 

0.313% 


;  -g-!400 

5  I ’.200 
s  »  ’000 
5  800 
E  —  600 


1  c 

’  642.3C0C0 

4  28*'. 

8  120 

0.256"'. 

1  ^ 

1  643.30000 

3  S3"'. 

7  934 

0.1  70".'. 

1  2 

1  536.80000 

3.29"'. 

7.710 

0.155% 

data.run«8 

t  (A.-gsiorm  ) 

S;c.Oe,._; 

mo  (10'ji_cm) 

Sl(J.Dev._fho 

1 

3500 

1 

2 

2860.80000 

1  0  08". 

7  355 

1  166% 

3300 

3 

2828  30000 

10  91  "i 

7.1  90 

1.2’5% 

e  3100 

4 

2473  80000 

9  03"'. 

7  474 

0  891% 

S  ~!900 

5 

2369  80000 

9.58"'. 

7.515 

1.073"'. 

i  f!700 

(J  B 

6 

2255  30000 

9  •  C"i 

7  6*3 

1.007"'. 

2  « :500 

' 

2203.30000 

8  82", 

7  667 

0.837% 

*■  c  ;3oo 

3 

2247  30000 

9  09". 

7  843 

1 .054% 

1  i’lOO 

9 

2250  30000 

7.77". 

7.853 

0  715"'. 

^  900 

1  0 

2243  'COCO 

6  0  8"'. 

7.913 

0.474% 

1 7C0 

*  * 

2253  30000 

6  6  0"'. 

8  003 

0  548". 

1  2 

2054  00000 

6  88"'. 

8.025 

0.595% 

1  SvO 

2  4  6  8 

Water  Number 


Run  #8 


4  6  8  1  0  1  2 

Water  Number 


()ata.run«9 


1  (Argsior-i) 

Std.Cov  _i 

rho  (10'Ji_cm) 

Sld.Dov._rho 

1 

2 

78*5  60000 

4  08"'. 

1  7.444 

0  1 08". 

3 

757  1 . 1  0000 

2  3  0"'. 

18  886 

0.032"; 

7800 

4 

742c  20000 

2  34"', 

1  9  678 

0  04  7?; 

„  7400 

5 

6i  '  3.00000 

2  45"'. 

1  8  872 

0  04  a".-; 

ff> 

6 

5975  10000 

3  03"'. 

1  6  581 

0  046?; 

1  ?  ■°°° 

7 

5291.70000 

4  02”'. 

1  4  753 

0  072". 

i  5  1600 

8 

5099  80000 

3  4  6"'. 

1  4  642 

0  063". 

9 

5298  10000 

2.28"', 

1  5  863 

0.038". 

^  r>i200 

c 

1  0 

5590  70000 

1.44”', 

1  7  532 

0.027?'. 

1  i.BOO 

1  1 

5762.60000 

1.02"'. 

1  8  942 

0  019?'. 

1  2 

521  7.00000 

0.93".i 

18  912 

0.014"/. 

6400 

Run  #9 


u 

i  •  J  "  '  I  '■ 

4  6  8 

Water  Number 


— r20  0 
:  19.5  E 
—  r  19  0  u 

-  18  5  J 

-  ’8.0, 
:17.5'o 
-  17.0  C 
:16  5  ^ 
:i6o  = 

I  -155  = 
\  -  IS  0  J 
Ij- 14.5  ! 

-4i4.0* 

12 


data.run*10 


Run  #10 


t  (Argstorm  ) 

Sta.Dev._1 

rho  (1Cr52_<im) 

Sta.Dev._rho 

1 

2 

3599  60C00 

1.89% 

9  553 

0.010% 

3 

3676,60000 

1  47% 

8  975 

0.0107'. 

4 

3742.00000 

1.697. 

8  732 

0.0067'. 

5 

3823.80000 

1.617'. 

8  813 

0.022% 

6 

3887.60000 

1.32% 

8  805 

0.0167'. 

7 

4032.00000 

1.17"'. 

8  834 

0.018"'. 

8 

4200.50000 

0.88"'. 

8  817 

0.0127'. 

9 

4367. 1  0000 

0.817'. 

8  760 

0.0127'. 

<■7 

V 

4531.50000 

0.91"/. 

8.755 

0.0147'. 

1 

44*6.30000 

1 .03% 

8  850 

0.0167'. 

z 

3855.50000 

1.627'. 

9  018 

0.024"/. 

I  (Arcsto'm  )  SlC  Dev,_t 


1 

3235  3CCC0 

7.3S"i 

2 

2955  'OOCC 

4,73"'. 

3 

2 . .560300 

4  1  2", 

4 

2708  30CC0 

3,38"i 

s 

2681.80000 

4  1  3"i 

e 

2638  30000 

4  08"', 

7 

262 '. 30000 

3  83"i 

8 

2704  50000 

3.30"', 

9 

28' 4.50000 

2  63"', 

0 

252*. 80000 

2  6  9"'. 

2925. 1  COCO 

3  24"', 

'  2 

2512.80000 

3.51"', 

1  (Argstorm  ) 

Sid.  Oev._ 

6'57  20000 

3  067', 

2 

601  2.20003 

2.79"'. 

3 

5283  63000 

2, 60"', 

4 

46'5 .70000 

3  2'". 

5 

4018  '  0000 

3.3  9"'. 

6 

3571.80000 

3  36"'. 

7 

3285,60000 

3  61". 

8 

31  77.60000 

3  38"', 

9 

3' 97,50000 

3,04"'. 

1  0 

3274.10000 

3,06°'. 

1  1 

3325  10000 

2,72"'. 

1  2 

31  03.30000 

4,41  % 

dala.runai  1 


rno  {10'^2_cm)  Sia.Dev._rho 


Run  #11 


6  830 
7.346 
7. 508 
7.440 

7  351 
7  465 
7.439 
7  519 
7.487 
7.479 
7.406 
7.232 


0.7317'. 
0.4S7% 
0.342% 
0.250% 
0  2857'. 
0.3097'. 
0,283% 
0.225% 
0.1717'. 
0.1947', 
0.199% 
0.228% 


Wafer  Number 


aata.ruri#12 


rho  (10'i_cm)  Sia.0ov._rho 


1  0  498 
10  822 
1  0  4C3 
9.739 
9  125 
8  819 

8  644 
8  649 
8  691 
8  703 
8  711 
8  255 


0,0837'. 
0.1297'. 
0.1047'. 
0.1417', 
0.195% 
0.197% 
0.222% 
0.2047', 
0  183% 
0.181% 
0,1637'. 
0.286% 


Run  #12 


n.'ita  from  Runs  10.  11.  and  12 


5  5 


Table  A.U10)-(121: 


daia.runaia 


Run  #13 


t  {AngstOrmj 

1 

2  2290.600CO 

3  2219.80000 

4  2163.10000 

5  2103.60000 

6  2029.10000 

7  1  975.10000 

8  1996.40000 

9  2072.60000 

I  0  2*  59.80000 

II  21 SS  1  0000 
1  2  1812.30000 


Sld.0av._i  rtio  (ic' 


6.83% 

5.27% 

5.12% 

5  51% 
6.30"; 

6  4  1*'. 
6.20% 
5.67% 
4.91% 
4.40«'. 
S.iei-i 


l_cm) 

Std.Dev.__rho 

7.570 

0.390% 

6.968 

0.386% 

7.1 11 

0.463% 

7  264 

0.564% 

7.331 

0.607% 

7.385 

0.575% 

7.341 

0.409% 

7  486 

0.460% 

7.469 

0.406% 

7.412 

0.268% 

7.610 

0.499% 

data.run«l4 


1  (Angsiorm  ) 

Std.Dev,_t 

rtio  (10'^_cm) 

Sld.Dew._fho 

1  6769.10000 

3.13-. 

8  698 

0. 105% 

2  7*29.80000 

4,74«i 

6.887 

0. 1  S27i 

3  653030000 

4  4  0% 

6.191 

C.228% 

4  5454.10000 

4,94% 

6  185 

0.374% 

5  4534.70000 

5.39% 

6  725 

0.379% 

6  3867  80000 

5  80-. 

6,923 

0.453% 

7  3470,60000 

6,13% 

7.236 

0.409*'. 

8  3436,30000 

5  58% 

7.471 

0  48B*i 

9  3548  60000 

5.03% 

7.519 

0  4267', 

’  0  36'8. 1  0000 

4,99% 

7. 507 

0.409% 

1  1  3648,  *  0000 

4.4S-, 

7  446 

0.290% 

12  3*S0, 30000 

5.33% 

6.895 

0.3247; 

daia.runal  s 

1  (Angsiorm  ) 

Std.Dev,_t 

rho  (1oi_cm) 

Std.Dev._rho 

1  6176.20000 

5  38% 

8.400 

0  4067; 

2  6885  70000 

5.72-. 

9  275 

0  3557; 

3  65«7  OCCCO 

4  43% 

9  519 

0.2837; 

4  6323  80CC0 

4,86% 

9  651 

0  215% 

5  6205.00000 

4,85% 

10.164 

0.2827'. 

6  6134,00000 

4.1  1% 

10.722 

0  230% 

7  6213.80000 

3  84% 

1  1.564 

0218% 

8  6398.20000 

3.37% 

12.310 

0.1737'. 

9  6647.30000 

3  25% 

13.168 

0.147% 

10  6707.20000 

2  05% 

13.927 

0.1207', 

11  6737,00000 

2.30% 

14  451 

0  0987'. 

12  5701.30000 

2.73% 

13.535 

0.103% 

Run  #14 


Wafer  Number 


Run  #15 


Wafer  Number 


Table  A,l(13)-(15):  Data  from  Runs  13.  14.  and  15 


56 


Reel»flvlly(10  J1  _cm|  ReelslWUy  _cm)  RealallvHy  (»0-‘'n_em) 


Run  #16 


Saia.run<l6 


t  (Angslirm  ) 

Sia.Dav._t 

rho  (iD'a-em) 

Sld.Dew._mo 

1 

461  7.00CCO 

3  07-. 

9,049 

0.054% 

2 

5257.10000 

3.03% 

8  522 

0,073% 

3 

5332.70000 

2.58r. 

8  591 

0,061% 

4 

5234.80000 

2.49% 

8  538 

0.033% 

5 

51  65.20000 

2,77% 

8.585 

0.04  2% 

6 

4980.00000 

2,39% 

8.322 

0.039% 

7 

4853.50000 

3  55-. 

8  023 

0,052% 

8 

4925,20000 

2.91% 

7.936 

0.040% 

9 

51  15.10000 

2  31% 

8.123 

0.027% 

1  0 

5254,70000 

1.7’% 

8  423 

0.014-. 

1  1 

52' 7  1 0000 

1 .39-i 

8  786 

0.027% 

1  2 

4483.20000 

2.25% 

8.307 

0.044% 

data.runflT 


t  (Angstcrn) 

Std,Dev._i 

rho  (10il_cfh) 

Sld.Dev._rho 

i 

4907  20000 

9  3 1 

5  864 

1.081% 

5059  80000 

9  98% 

S.9C5 

1.421% 

3 

433' ,00000 

7.07-. 

6.263 

0.949% 

4 

3891  80000 

6  83% 

6.674 

0  644% 

5 

3499  3CC00 

6  82-. 

7.072 

0.658% 

6 

3229  30000 

7.31V. 

7.180 

0.805% 

7 

2383  60000 

8.03-. 

7.507 

1.013% 

3 

2952  50000 

8.05% 

7.538 

1.000% 

9 

2379  00000 

6,71-. 

7,471 

0.703% 

1  0 

3036  30000 

6,46-. 

7.339 

0  439% 

1  1 

3003  80000 

6  18-. 

7,488 

0.590% 

1  2 

2559.30000 

6.93-. 

7.643 

0.716% 

Run  #17 


aala.run«18 


t  (Angsicrm)  Sia  Cev.J  rtio  (10jl_cmj  Std.Dev._rho 


1  5251. 

2  5:56. 
53<i: 
50'6 
4876 
4746 

7  4632, 

8  47  55 

9  4822 

10  4846 
1  1  4830 

12  4190 


70C00 
70CC0 
3CC00 
20000 
SCCCO 
1  0000 
20000 
20000 
.70000 
.80000 
20000 
.00000 


8  0 1  r. 

8.92-. 
8  40% 
8.36-. 
8.50% 
8  61-. 

8  47% 

9  59% 
7.55% 
7,02% 
6.85% 
7.25% 


6.843 
7  035 
6  914 
6,907 

6  915 
7.005 
7.171 
7.442 

7  321 
7.333 
7.385 
7,253 


0  812% 

1 .237% 
1.11 0% 
1.172% 
1.153% 
1.178% 
1.017% 
1.401% 
0.866% 
0.896% 
0.709% 
0.753% 


Run  #18 


Wafer  Number 


Table  A.HI6)-(18):  Data  from  Runs,] 6,  17.  and  JJ 
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cm>  nn*  fl  rmy 


Degrees  of 
Factor  Freedom 


Sum  of 
Squares 


Meart  Sum  of 
Squares 


% 

Effect 


Rank 


F-Ratio 


35.422 


176.405 


72.430 


14.944 


29.563 


5.619 


309.799 


a2 

2 

149.514 

Error 

2 

1.793 

Total 

1  7 

795.487 

35.422 


88.202 


36.215 


7.472 


14.782 


2.809 


154.899 


74.757 


0.897 


39.508 


98.377 


40.393 


8.334 


16.487 


3.133 


72.768 


83.381 


4.5% 


22.2% 


9.1% 


1 .9% 


3.7% 


0.7% 


18.8% 


Degrees  of 

Sum  of 

Mean  Sum  of 

Factor 

Freedom 

Squares 

Squares 

W 

1 

.789 

.789 

HB 

2 

38.399 

19.199 

T 

2 

8.092 

4.046 

P 

2 

2.498 

1.249 

■B 

2 

35.763 

17.882 

Qs 

2 

22.296 

11.148 

Qd  f 

2 

39.240 

19.620 

ES 

2 

219.856 

109.928 

Error 

2 

43.391 

21.696 

Total 

1  7 

410.323 

a2 


Error 


Total 


Degrees  of 

Sum  of 

Freedom 

Squares 

1 

6.040 

2 

1.151 

2 

13.553 

2 

0.068 

2 

3.189 

2 

2.886 

2 

32.609 

2 

15.343 

2 

0.340 

1  7 

75.178 

Mean  Sum  of 
Squares 


6.040 


0.576 


6.777 


0.034 


1.595 


1.443 


16.304 


7.671 


0.170 


F-Ratio 


35.483 


6.761 


79.625 


0.398 


18.737 


16.953 


191.577 


90.139 


% 

Effect 


8.0% 


1 .5% 


18.0% 


4.2% 


3.8% 


43.4% 


20.4% 


Table  A. 4:  ANOVAJor  SNq 
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Degrees  of 

Sum  of 

Mean  Sum  of 

Factor 

Freedom 

Squares 

Squares 

W 

1 

38.334 

38.334 

Ac 

2 

20.305 

10,152 

T 

2 

5.682 

2.841 

P 

2 

10.742 

5.371 

mm 

2 

27.982 

13.991 

Qs 

2 

11.352 

5.676 

Qd1 

2 

61.881 

30.941 

E9 

2 

32.799 

16.400 

99SB 

2 

12.739 

6.369 

iTotal 

1  7 

221.816 

F-Ratio 


6.019 


1.594 


0.446 


0.843 


2.197 


0.891 


4.858 


2,575 


% 

Effect 


17.3% 


9.2% 


2.6% 


4.8% 


1 2.6% 


5.1% 


27.9% 


14.8% 


t  (Angslorm)  Sid.  Dev._l 


1  4952  2000C0 

2  4440  500C00 

3  4396.100000 

4 

5  4463  000000 

6 

7  4589.000000 
6 

9  4877.300000 
1  0 
1  1 

1  2  3983  500000 


0  787V. 
0  954% 

0  334% 


1  (Angslorm)  Sid.  Dev._l 


3  2921.300000 

4 


9  31  72.100000  1.39B5-. 

1  0 

1  1  3254  100000  0  733". 

1  2  3855.100000  0  5ei5; 


I  (Angslorm)  Sid  Osv  _( 


2  5139.500000  4  66f.i 

3 


8  5703.500000  1.363*'. 

7 

8 

9 

1  0 

1  1  651  8.700000  0  OlO’.i 


rho  (idlt  i 


1 

2  12  330669 

3  12.710015 

4 

5  12881547 

8 

7  13.070078 

8 

9  12  577869 

1  0 

11  12  239948 

1  2  13.490182 
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J  o  2500 
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There  are  two  kinds  of  film  thickness  quality  loss  functions,  across- 
the-wafer  thickness  variance  cost,  and  the  cost  due  to  the  total 
variance  across  a  batch  (across-the-wafer  "plus"  down-the-load): 

(a)  Quality  Cost  due  to  Across-the-wafer  Variance: 

Each  batch  consists  of  nw  wafers  with  and  "Ox,i"  as  the  mean 

and  standard  deviation  of  the  film  thickness  of  the  ith  wafer.  Quality 
loss  due  to  ith  wafer  in  a  batch  is  the  sum  of  the  squared  deviation  of 
the  average  value  of  wafer  film  thickness  from  target  and  the 

resulting  mean  squared  deviation  of  film  thickness  around  its  own 
mean. 


(B.l) 


where  kxi  is  the  corresponding  quality  loss  coefficient  that  converts 
the  loss  into  a  S  value.  Now'  if  a  adjustment  factor  like  deposition 
time  is  increased  by  a  factor  (^,i/^T,i)  so  as  to  obtain  the  target 
thickness,  the  corresponding  standard  deviation  after  adjusting  the 
mean  changes  to  and  the  adjusted  quality  loss  for  a 

wafer  is  given  by 


'^T.i 


(B.2) 


The  total  adjusted  quality  loss  for  a  batch,  as  defined  by  Taguchi,  can 
be  expressed  in  two  ways.  Firstly,  the  sum  of  losses  incurred  by 
individual  wafers  (Equation  B.2)  can  be  used  to  represent  the  total 
average  loss  for  a  batch.  Note  that  this  loss  is  only  due  to  the 
thickness  variance  across  the  wafers  and  is  given  by 


riw 


(B.3) 


In  Equation  B.3  it  is  assumed  that  the  quality  loss  coefficient  kji  is 
same  for  all  wafers  in  a  batch.  This  is  justified  because  in  case  a 
batch  has  to  be  scrapped,  all  the  wafers  have  same  dollar  expense 
inccured  in  them. 


(b)  Quality  Cost  due  to  Total  Batch  Variance: 

Once  again,  considering  that  each  batch  consists  of  nw  wafers  with 
and  as  the  mean  and  standard  deviation  of  the  film 

thickness  of  the  ith  wafer,  the  overall  mean  for  the  batch  is  given  by 


'^1=1 


(B.4) 


If  the  thickness  variance  of  the  ith  wafer  is  based  on  m,  data  points 
across  a  wafer  then  the  batch  variance  is  given  by 


S  +4.,K 


<4=— 


Im, 

i=l 


(B.5) 


In  equation  B.5  the  variance  across  the  wafer  Oj  j  is  "added"  to  the 
variance  of  the  mean  l^T.i  down  the  load.  If  the  variance  of  all  wafers 
is  based  on  same  number  of  data  points  then 


isl 


(B.6) 


Therefore,  the  total  batch  variance  incorporating  down-the-load 
variance  as  well  as  across-the-wafer. variance  is  given  by. 


62 


(B.7) 


Q^T:  ~  ^T2 


where  kj:  is  the  loss  coefficient  for  total  variance  in  a  batch. 

For  a  target  film  thickness  of  ^r.t  (A°)  consider  Aq  as  the  S  cost 
incurred  in  discarding  a  batch  if  across-the-wafer  nonuniformity 
exceeds  a  functional  limit  of  "Vi  (A°).  Using  this  in  Equation  B.4.  we 
get  , 


(B.8) 


Similarly,  imposing  the  functional  limit  of  ‘^2  (A°)  on  total  thickness 
variation  in  batch  (across-the-wafer  "plus"  down-the-load).  we  get. 


(B.9) 


(2)  Growth  Rate  Loss: 


Growth  rate  degradation  is  not  a  quality  loss  because  to  compensate 
for  lower  growth  rate  one  can  increase  the  deposition  time  without 
scrapping  the  lot.  Therefore,  the  loss  incurred  due  to  increased 
deposition  time  is  only  the  increased  total  labor,  equipment,  and 
floor  space  costs  Cl,  Ce,  Cf  (S)  respectively.  For  a  w-afer,  this  cost  falls 
as  the  average  growth  rate  for  a  wafer  rises.  Again,  summing  the 
individual  wafer  costs,  we  get  the  average  batch  cost  for  lower 
growth  rates  as. 


Cc 


~  Pt  1 


(B.IO) 
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where  ko  is  the  loss  coefficient  that  expresses  the  above  loss  in  S 
value.  Consider  a  fixed  cost  rate  cl.ce.  and  cp  ($/min)  for  labor, 
equipment,  and  floor  space  respectively.  The  total  cost  incurred  for 
deposition  time  to  is 

CG=(to+tl)(CL+CE+CF)  (B.ll) 


where  ti  is  the  fixed  time  for  predeposition  and  postdeposition  steps. 
If  the  deposition  time  to  is  increased  to  then  the  film  thickness 
achieves  the  target  Using  this  fact  in  Equation  B.ll  and  after 

readjustment,  we  get. 


k^  =(CL+Cp+q) 


r  *» 

Lmin-' 


(B.12) 


(3)  Film  Resistivity  Quality  Loss; 


The  target  is  to  minimize  the  resistivity  of  the  film.  For  zero  quality 
loss  the  desired  resistivity  is  zero.  However,  physically  it  is 
impossible  to  obtain  zero  resistivity.  Therefore,  any  value  of 
resistivity  leads  to  some  loss.  Modelling  this  loss  as  a  quadratic  loss, 
we  get  the  average  loss  for  a  batch  as 


(B.13) 


where  kp  is  the  resistivity  quadratic  loss  coefficient.  Now'  consider 
that  the  batch  is  scrapped  at  a  loss  of  Aq  ($)  if  the  average  film 

resistivity  in  a  batch  exceeds  the  functional  limit  of  above  the 
target  of  Fp  Q_cm  (which  is  zero  ideally).  Using  these  facts  in 
equation  B.13,  we  get. 


k 


{B.14; 
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